Spin Anisotropy Effects in Dimer Single Molecule Magnets 
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We present a model of equal spin si dimer single molecule magnets. The spins within each dimer 
interact via the Heisenberg and the most general set of four quadratic anisotropic spin interactions 
with respective strengths J and {Jj}, and with the magnetic induction B. We solve the model 
exactly for si = 1/2, 1, and 5/2, and for antiferromagnetic Heisenberg couplings (J < 0), present 
M(B) curves at low T for these cases. Low-T Cv(B) curves for si = 1/2 and electron paramagnetic 
susceptibility \(B,lo) for s\ = 1 are also provided. For weakly anisotropic dimers, the Hartree 
approximation yields rather simple analytic formulas for M(B) and C'v(B) at arbitrary Si that 
accurately fit the exact solutions at sufficiently low T or large B. Low-T, large-B formulas for 
the inelastic neutron scattering cross-section S{B,q,ui) and x(B,uj) with arbitrary si and B in an 
extended Hartree approximation are also given. We show that the various anisotropy interactions 
give rise to new, detectable electron paramagnetic resonances and to new peaks in S(q, u), which can 
provide precise determinations of the particular set of microscopic anisotropy interaction strengths. 
For antiferromagnetic Heisenberg couplings (J < 0) and weak anisotropy interactions (| Jj/J\ <C 1), 
we provide analytic formulas for the 2si level-crossing magnetic inductions Bf s (6, (f>), at which 
the low-T magnetization M(B) exhibits steps and the low-T specific heat Cv(B) exhibits zeroes, 
surrounded by double peaks of uniform height. Strong anisotropy interactions drastically alter these 
behaviors, however. Our results are discussed with regard to existing experiments on si = 5/2 Fe2 
dimers, suggesting the presence of single-ion anisotropy in one one them, but apparently without any 
sizeable global spin anisotropy. Further experiments on single crystals of these and some si = 9/2 
[Mn4]2 dimers are therefore warranted, and we particularly urge electron paramagnetic resonance 
and inelastic neutron scattering experiments to be performed. 

PACS numbers: 05.20.-y, 75.10.Hk, 75.75.+a, 05.45.-a 



I. INTRODUCTION 



Single molecule magnets (SMM's) have been under in- 
tense study recently, due to their potential uses in mag- 
netic storage and quantum computing. 0, 0, The ma- 
terials consist of insulating crystalline arrays of identical 
SMM's 1-3 nm in size, each containing two or more mag- 
netic ions. Since the magnetic ions in each SMM are 
surrounded by non-magnetic ligands, the intermolecular 
magnetic interactions are usually negligible. Although 
the most commonly studied SMM's are the high-spin 
Mni2 and Fes , III- L3- lil Lit LJ such SMM's contain a vari- 
ety of ferromagnetic (FM) and antiferromagnetic (AFM) 
intramolecular interactions, rendering unique fits to a va- 
riety of experiments difficult. 6] 

In addition, there have been many studies of AFM Fe„ , 
ring compounds, where n = 6, 8, 10, 12, etc.Q, |& Oj, lljj 
In these studies, analyses of inelastic neutron diffrac- 
tion data and the magnetic induction B dependence of 
the low-temperature T specific heat Cy and magneti- 
zation M steps were made, using the isotropic Heisen- 
berg near-neighbor exchange interaction, the Zeeman 
interaction, and various near-neighbor spin anisotropy 
interactions. |t1 liL 1^. Ho| However, the rings were so com- 
plicated that analyses of the data using those simple 
models were inaccessible to present day computers. 0, 
Thus, those authors used either simulations or phe- 



nomenological fits to a first-order perturbation expansion 
with different spin anisotropy values for each global ring 
spin value. H ll E3 

Here we focus on the simpler cases of equal spin s± = S2 
magnetic dimers, for which the full spin anisotropy effects 
can be evaluated analytically, investigated in detail nu- 
merically, and compared with experiment. AFM dimers 
with si = 1/2, 3/2,[I3 HG1 

LLJI and various forms of 
Fe? with si = 5/2 were studied recently. [l5L Hr j Il7l liH 
[ia|20( Several Fe2 dimers and effective si = 9/2 dimers 
of the type [Mru]9.|2ll |22j have magnetic interactions 
weak enough that their effects can be probed at T w IK 
with presently available B. A comparison of our results 
with M(B) step data on a Fe2 dimer strongly suggests 
a substantial presence of local spin anisotropy. ^3 

The paper is organized as follows. In Section II, we 
present the model in the crystal representation and give 
exact formulas for the matrix elements. The general ther- 
modynamics are presented in Section III, along with the 
universal behavior of the M(B) and Cy{B) behavior as- 
sociated with the energy level crossing. In Section IV, 
we solve the model exactly for s\ — 1/2, giving analytic 
expressions for M and the specific heat Cy. In Sec- 
tion V, we discuss the exact solution for s\ = 1, present 
the equations from which the eigenvalues are readily ob- 
tained, and give examples of the low-T M(B) curves. In 
Section VI, we present examples of the low-T si = 5/2 
M{B) curves. In Section VII, we present our induction 
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representation results for the eigenstates to first order in 
the anisotropy energies. These are used to obtain the 
asymptotic Hartree expressions for the thermodynamic 
quantities M and Cy , and for the inelastic neutron scat- 
tering cross-section S(B,q,u) and the electron param- 
agnetic resonance susceptibility x(-B, w) for arbitrary s±, 
that are accurate at sufficiently low T and/or large B. 
In addition, an analytic formula accurate to second or- 
der in the four independent anisotropy exchange energies 
is provided for B' C S1 (0, </>), the B values at which the sth 
energy level crossing occurs. In Section VIII, we consider 
FM anisotropic interactions sufficiently strong as to re- 
move the level crossing effects. Finally, in Section IX, we 
summarize our results and discuss them with regard to 
experiments on Fe2 dimers. 



II. THE MODEL IN THE CRYSTAL 
REPRESENTATION 

We represent the s± = S2 dimer quantum states, \ipT) 
in terms of the global (total) spin and magnetic quantum 
numbers s and m, where S = Si + S2 and S z = S ■ z 
satisfy S 2 \ip™) = s(s + 1)|V>™) and Sz\4>?) = "#D> 
where s = 0, 1, . . . , 2s\, m = — s, . . . , s, and we set h = 1. 
We also have S±\ip?) = Af" 1 ^^ 1 }, where S± = S x ± 
iS y and 



a: 



y (s — m)(s + m + 1). 



(1) 



For an arbitrary B, we assume the Hamiltonian has the 
form 7i = Hq + H a + Kb + Hd + H e , where 



Ho = -JS 2 /2-jS-B 



(2) 



contains the Heisenberg exchange and Zeeman interac- 
tions, the gyromagnetic ratio 7 = g^B, 9 ~ 2 and \ib 
is the Bohr magneton. The global axial and azimuthal 
anisotropy terms 



Hb 



-JbS, 



and 



Hd — Jd{S x 



(3) 



(4) 



respectively, only involve components of S, but have 
been the main anisotropy terms discussed in the SMM 
literature, 0, 0| so we have included them for compari- 
son. Such terms have been commonly studied as an ef- 
fective Hamiltonian for a singlet orbital ground state, in 
which the tensor global spin interaction with a fixed spin 
quantum number s has the form S ■ A • S, resulting in the 
principal axes x,y and z.|2.lj For a dimer, we take A to 
be diagonal in the orientation pictured in Fig. 1. Then 
J = -{A xx + Ayy)/2, J b = -A 22 + (A xx + Aj,„)/2, and 
Jd = (A yy - A xx )/2. Taking \ J d /J\ < 1 and \ J b /J\ < 1 



still leaves Jd/Jb unrestricted. The single- ion axial and 
azimuthal anisotropy terms, 



and 



H a J a ^ ' "5*i2 
i=l 

2 

= —Je ^ (Sj x — Sjy 



t=l 



(5) 
(6) 



respectively, arise from spin-orbit interactions of the local 
crystal field with the individual spins. These terms have 
usually been neglected in the SMM literature, but have 
been studied with regard to complexes containing a single 
magnetic ion, such as Ni +2 ,|2^| and with regard to clus- 
ters of larger numbers of identical magnetic ions.[23ll2^| 
The local axially and azimuthally anisotropic exchange 
interactions 



Hf — —JfSi z S2 Z , 

H c = —J c (£>lxS2x ~ SlyS2y), 



satisfy 



Wf/Jf 
2H C /J C 



Hb/ Jb — Hal Ja, 
Hd/ Jd — H e / Je, 



(7) 
(8) 



(9) 
(10) 



so we need only include either either H a or Hf and H e 
or H c , respectively. [23| That is, if we stick to the Hamil- 
tonian H, we may incorporate the effects of Hf and 
H c by letting J b — > J b + J//2, J a -> J a - J//2, and 
Jd — » Jd + J e/ 2, Je — > J e - Jc/2, respectively. Since H a 
and He describe the axial and azimuthal anisotropy each 
single ion attains from its surrounding environment, they 
are the physically relevant local anisotropy interactions. 
We note that Hb and H a are symmetric under x <-> y, 
whereas Hd and H e are antisymmetric under x <-> y, 
independent of s±. 

For the case of Fe2,^3 a constituent of the high-spin 
SMM Fe 8 and the AFM Fe n rings,[lH0] the exchange 
between the Fe +3 si = 5/2 spins occurs via two oxygen 
ions, and these four ions essentially lie in the same (xz) 
plane. [l^Lllsj l We set the z axis parallel to the dimer axis, 
as pictured in Fig. 1. Since the quantization axis is along 
the dimer axis, which is fixed in a crystal, we denote this 
representation as the crystal representation. 

We generally expect each of the Jj for j = a, b, d, e to 
satisfy \Ji/J\ <C 1, but there are not generally any other 
restrictions upon the various magnitudes of the Jj . Since 
all dimers known to date have predominantly AFM cou- 
plings ( J < 0), and also because their magnetizations and 
specific heats are particularly interesting, we shall only 
consider AFM dimers. In addition, since no studies on 
unequal-spin dimers have been reported to our knowl- 
edge, we shall only treat the equal-spin si = S2 case. 
We note that for equal spin dimers, the group symmetry 
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and 



FIG. 1: Sketch of an Fe2 dimer, with two bridging -2 ions 
(O). Ligands (not pictured) are attached to the Fe +3 ions 
(Fe). The arrows signify spins. 



of the dimer environment is Ci v , so that Dzaloshinskii- 
Moriya interactions do not arise, jl^ Hence, our Hamil- 
tonian H is the most general quadratic anisotropic spin 
Hamiltonian of an equal-spin dimer. 

For B = B(sin 9 cos <f>, sin 9 sin <fi, cos 9), we have 



n \r s 
n b \r s 

and 

n d \r s 

where 



<T = ±1 



E™\if)™)+8E e- Ur +A'? n \il>? +c 

-VIC), 



cr=±l 



E? 
SE 
b 

and 

para 



— Js(s + l)/2 — mbcost 

— 6sin#, 
2 

jB. 



A am a l+erra 



(11) 

(12) 



(13) 



(14) 

(15) 
(16) 

(17) 



TL a and H. c contain the individual spin operators Si Z 
and Si± for i = 1,2. After some Clebsch- Gordon alge- 
bra involving the Wigner-Eckart theorem, for arbitrary 
(s 1 ,s,m), 



Si±W s 



l_ A ± mw±1) T I(_i)i(ci£|^±i) 

-cjIITM?^)), (is) 
+dt +1 ,sMT + i)), (i9) 



where 



Dl 



Sl y/(s- m)(s -m- 1), 
si V / ( fi2 - ™ 2 ), 



(20) 
(21) 



^[( 2si + l)2_ s 2 ]/(4s 2_ 1) . 



(22) 



The generalizations of these results to dimers with gen- 
eral s\, S2 are given in Appendix A. 

For s = 0, we require m = 0, for which Cq s = Dp si = 
0. We then find 



er'=±l 

<T=±1 

+ £ ^T; ff 'lC + + 2 2 CT ' 



Je 

4 



G 



and 



S,Sl 



Tjiyi.a 
H s,si 



<T>=±1 



s(s + 1) - 1 

+ [2m 2 + l-2s(s + l)]a SjSl , 

ly s+(a' + l)/2,s 1 ly s+(3cr' + l')/2,s 1 > 

r ,-<T'a;-(l+ ( T')/2 r ,-o-'2;-(3cr' + l)/2 
Ly s+(<r'+l)/2,si U s+(3it'+1)/2,si ' 

3s(s + 1) - 4si(si + 1) - 3 
(2s- 1) (2s + 3) ' 



(23) 



(24) 



(25) 
(26) 

(27) 
(28) 

(29) 



We note that 2a S}S1 = 1 — rjl — ?7j +1 a and that Eq. 
O holds for s > 0. Equations JTBJ) and JTSJl allow 
for an exact solution to the most general Hamiltonian of 
arbitrary order in the individual spin operators. In all 
previous treatments of more complicated spin systems 
with similar anisotropic interactions, it was only possi- 
ble to obtain numerical solutions, and therefore the full 
anisotropy of the mag netization and specific heat was 
not calculated. [i^, l24| The operations of Ho, Hb and 
Hd satisfy the selection rules As = 0, Am = 0, ±1,±2. 
The local anisotropy interactions H. a and H. e allow tran- 
sitions satisfying As = 0,±2, Am = 0, and As = 0,±2, 
Am = ±2, respectively, so in the presence of either of 
these interactions, s is no longer a good quantum num- 
ber, unless si = 1/2. 



III. GENERAL THERMODYNAMICS 

In order to obtain the thermodynamic properties, we 
first calculate the canonical partition function, Z = 
Tiexp(—f3H). Since H. is not diagonal in the (s,m) rep- 
resentation, we must construct the wave function from 
all possible spin states. We then write 



(30) 
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where \^ Sl ) is constructed from the {IV'™)} basis as 

<* s j = ((^i.<^."r 1 i.-.^?i.(c 1 i.^oi), 

(31) 

where /? = l/(kgT) and kg is Boltzmann's constant. 
To evaluate the trace, it is useful to diagonalize the 
(* ai |ft|tf ai ) matrix. To do so, we let |* Sl ) = U\$ ai ), 
where 



= ((0n sl l,(</V n -l|,---,(</>l|): 



(32) 



is constructed from the new orthonormal basis {|<^ n )}, 
and U is a unitary matrix of rank n si — (2si + l) 2 . 
Choosing U to diagonalize H, UHW = H, we generally 
obtain 7i\(f> n ) — e n \<j) n ) and the partition function for a 
SMM dimer, 



Z = ^exp(-/3e„) 



(33) 



n=l 



The specific heat Cy = kg^d 2 \nZ/d{3 2 is then easily 
found at all T, B, 



e 2 e'^ 



(E 



-/3e n 



.(34) 



Z 2 

The magnetization 

M = --^V B (e„)exp(- / 3e„), (35) 

71=1 

requires Vs(e n ) for each B. As T — > 0, at most two 
eigenstates are relevant. For most B values, only one e n 
is important. But near the reduced sth level-crossing in- 
duction 6* = r yB\ c Si (9,(j)) at which e s = e s _x, two eigen- 
states are relevant. We then obtain for these two energies 



C v (b)/k B 
M(b) 



(36) 



(/3Ae s /2) 2 
cosh 2 (/3Ae s /2)' 

i(M s + M s _ x ) 

+1(M S - M._i) tanh(/3Ae s /2), (37) 



where 



Ae, 



e.(6)-c._i(6), 



(38) 



and M s = 7V{,e s (&) is the magnetization of the sth eigen- 
state only. We then expand in powers of b — b* , 



e a (b) = e s (b:) + (b-b:) ai , s + ±(b-K) 2 a 2 , s + 



and we find 



c v (K) 



0, 



(39) 
(40) 



Cy 



b*± 



2c 



— dl,s-l 

C v cak /k B 



k B T<^\J\ 



c 



peak 



(41) 



(-Y 

Vcoshc/ 
0.439229, (42) 



and 



M(K)/j T Z -(a M +a M -i) 



dM 



= a -1/2 + 0(4/ J) 2 , (43) 
= |(l + 0(4/ J) 2 ), (44) 



where s = 1, . . . , 2si and c 1.19967864 is the solution 
to tanhc = 1/c. The coefficients di )S = s + 0(4/ J) 2 
and a2,s = 0(Jj/J) 2 . The easiest way to see this is to 
first rotate the crystal so the quantization axis is along 
B, as discussed in Appendix B. An expression for ai )S 
to second order in the Jj/ J is given in Appendix D. We 
note that b* = r yB 1 f Si (9,4>) depends upon s,si, and the 
direction of B when anisotropic interactions are present. 

The Cy(B) double peaks are double Schottky anoma- 
lies at the reduced induction values b* + + 0(Jj/J) 2 . 
Hence, to (Jj/J), the heights and midpoint slopes of the 
2si low-T M(B) steps are uniform and the same as for 
the isotropic case, but the step positions and hence their 
plateaus are not. Correspondingly, the heights and po- 
sitions of the 2si Cy(B) double peaks are uniform and 
the same as for the isotropic case, but positions of the 
zeroes are not. Hence, the non-universal level-crossing 
inductions B° (9, <ft) fully determine the low-T thermo- 
dynamics of weakly anisotropic AFM dimers. 

In the next three sections, we consider the special 
cases of Si = 1/2,1 and 5/2. Then, in Section VII 
and Appendix D, we present out general expression for 
B]f (9, 4>) accurate to second order in each of the Jj . We 
remark that a double peak in the low-T Cy(B) curve 
has been seen experimentally in a much more compli- 
cated Fen ring compound, and was attributed to level 
crossing, yj 



IV. ANALYTIC RESULTS FOR SPIN 1/2 

Plots of Cy/ks and M/7 versus jB/\J\ for the 
isotropic spin 1/2 dimer were given previously. |l4j| For 
s% = 1/2 with an arbitrary B and Jj for j = a, b, d, e, 
the rank 4 Hamiltonian matrix is block diagonal, since 
s = 0, 1 is a good quantum number. The eigenvalues are 
given by 



ei = - 



Ja 

2 ' 

Ja 

2 



J + A n , 



(45) 

2,3,4, (46) 
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where 







-b 2 sin 2 9[J b - J d cos(20)]. (47) 
For the special cases B\\i for i = x,y, z, the \ l n satisfy 



\ z n — 0,—Jt,±F z 

yx,y _ _ 2 7 



-J + F 



where 



Ft = \Jb 2 + Jl 

Jx,y,z = (Jb ± Jd)/2, Jd, 



(48) 
(49) 



(50) 
(51) 



respectively, and where J x (J y ) corresponds to the upper 
(lower) sign. 

The magnetization for B\\i is given by 

76sinh(/3Fi) , . 

M t = 7 ^ l \ (52) 

Vi = cosh((3Fi) + Ai/2, (53) 

A, = exp(-/3J :c )[exp(-/3J)+exp(2/3J y )], (54) 

A y = exp(-/3J y )[exp(-/3J)+exp(2/3J 2; )], (55) 



exp(-/3J 6 )[exp(-/3J) + 1], 



(56) 



respectively. When the interactions are written in terms 
of the less physical Jf,Jb,Jd, and J c , then Jb — > J b + Jf 
and Jd ^ Jd + Jc/2, so that for si — 1/2, Jf and J c 
merely renormalize Jb and Jd- Neither of the single- ion 
spin anisotropy terms H. a and H. e affect the thermody- 
namics for si = 1/2. We note that M y (J,i) = M x (— Jd) 
for each J5, as expected from H. d = —Jd{S 2 — S 2 ). 

From Eqs. (JEJl, (001, EEl, and igHJ), the single level 
crossing induction (s = 1) for am S\ = 1/2 dimer with 
B\\i occurs at 



1.0 



S=1/2 AFM 



to 



k B T/|J|=0.03 / / 



J =0.1 J 

b 




9M B B/|J| 



FIG. 2: Plots of M/7 versus jB/\J\ at k B T/\J\ = 0.03 for 
the AFM spin 1/2 dimer with J b = 0.1J, J d = 0, with B\\z 
(solid), B\\x (dashed), along with the isotropic case Jb = 
J d = (dotted). 




0.75 



1.00 



gn B B/|J| 



1.25 



1/2 



\J J 2 + J{Jb ± Jd) 



Jb) 2 



p 



B\\x,y 
B\\z 



(57) 



provided that J + J x y < and J + Jb < 0, respectively. 

To distinguish the different effects of the global 
anisotropy interactions Jb and J d that affect the mag- 
netization of s\ = 1/2 dimers, in Figs. 2 and 3, we 
have respectively plotted the low-T M/7 versus r )B/\J\ 
with J b = 0.1J, J d = and J d = 0.1 J, J b = for B|||z 
(solid) and B\\x (dashed), along with the isotropic case 
J b = Jd = (dotted). From Fig. 2, J b < and J d = 
causes a greater shift to higher B values at the magneti- 
zation step with B\\z than for .B||i, consistent with Eq. 
(|57|1 . In addition, J d finite with Jb = has a very differ- 
ent effect upon the anisotropy of the magnetization step, 
as shown in Fig. 3. Although for B\\z, B at the step is 
slightly reduced from its isotropic interaction value, for 
£?||a:, the magnetization step occurs at a larger B. These 



FIG. 3: Plots of M/7 versus jB/\J\ at k B T/\J\ = 0.03 for 
the AFM spin 1/2 dimer with J d = 0.1J, J b = 0, with B\\z 
(solid), B\\x (dashed), along with the isotropic case Jb = 
J d = (dotted). 



results are consistent with Eq. I|57l) . The midpoint slopes 
are universal, in accordance with Eq. 144fl . 

The specific heat of an si = 1/2 dimer with B\\i is 



Cvi — 



k B P 2 K 



(58) 



where the T>i are given by Eq. (|53|l . and the Mi are 
given in Appendix A. Plots at low T of Cv/ks versus 
r )B/\J\ for si = 1/2 dimers with the corresponding global 
anisotropics Jb — 0.1 J and Jd — 0.1J are shown in Figs. 
4 and 5, respectively. We note the universal curve shapes, 
but non-universal level-crossing positions, in quantitative 
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FIG. 4: Plot of C v /k B versus jB/\J\ for the AFM spin 1/2 
dimer with J h = 0.1J, J d = at k B T/\J\ = 0.03 with the 
same curve notation as in Fig. 2. 



O 0.2 



S=1/2AFM 




Jd=0 ' 1J A' / 


V\ 




1\ 1^1/^1=0.03 


j ' i J 


'\ 1 
i i i 


biiz^/ ; 1 i 

y ' 1 J 
/ * If 




y i i ] 

J / In 
/i Hi, 


\ I 
\ \ 
\ \ 
\\ 



9h b b/|J| 

FIG. 5: Plot of C v /k B versus jB/\J\ for the AFM spin 1/2 
dimer with J d = 0.1J, J b = at k B T/\J\ = 0.03 with the 
same curve notation as in Fig. 3. 



agreement with Eqs. H40() - (|42() . In Fig. 4, the positions 
of the maxima and the central minimum in Cy track 
that of the magnetization step in Fig. 2 with the same 
parameters. With Jd = 0.1J, the behaviors in Cv and 
M for B\\x are also very similar. However, there is a 
slight difference in the behaviors for B\\z. Note that M 
(Fig. 3) shows a slight reduction for B\ \z in the induction 
required for the step, whereas CV (Fig- 5) shows a slight 
increase in the positions of the peaks. This detail only 
appears when Jd ^ 0, for which the effective temperature 
is slightly higher than for Jd = 0. 



V. ANALYTIC AND NUMERICAL RESULTS 
FOR SPIN 1 

For dimers with S\ = 1, the allowed s values are s = 
0, 1, 2. The three s = 1 states are decoupled from the six 
remaining s = 0, 2 states. They satisfy a cubic equation 
given in Appendix A. For B\\i with i = x,y, z, this cubic 
equation simplifies to a linear and a quadratic equation, 
as for the s = 1 eigenstates of si = 1/2 dimers. 

The remaining six eigenstates corresponding nominally 
to s = 0, 2 are in general all mixed. The matrix leading to 
the hexatic equation from which the six eigenvalues can 
be obtained is given in Appendix A. That is sufficient 
to evaluate the eigenvalues for si = 1 using symbolic 
manipulation software. When combined with the three 
s = 1 eigenvalues, one can then use Eqs. (|34|l and l|35|) 
to obtain the resulting exact magnetization and specific 
heat at an arbitrary B. The combined nine eigenvalues 
depend upon all four anisotropy parameters Jj for j = 
a, b, d, e. 

To the extent that the eigenvalues can be obtained 
from the solutions to either linear or quadratic equations, 
the expressions for the B l f 1 are simple, and are given in 
Appendix A. In the global anisotropy case J a = J e = 0, 
the first level crossing induction B l f 1 with B\\i is iden- 
tical to B 1 ^^, the level crossing with s% = 1/2 given by 
Eq. I|57p. For comparison, we expand the first and sec- 
ond level crossing inductions to first order in each of the 
Jj for j — a, b, d, e for B\\i where i = x,y, z, 



7*S2 = 


-J + 


Ja 

3 


Jb, 




(59) 


7<:l = 


-J - 


Ja 

6 


Jb l fT 

Y T 2 {Jd 


" Je), 


(60) 


7<S = 


-2J- 




-3J 6 , 




(61) 


7<L = 


-2J- 




Jb 5J d 
2 T 2 


3J e 


(62) 



where the upper (lower) sign is for B\\x (B\\y), respec- 
tively. 

From these simple first-order results, it is possible to 
understand the qualitatively different behavior obtained 
with local, single-ion, anisotropy from that obtained with 
global anisotropy. In the isotropic case Jj — OVj, the first 
and second level crossings occur at — J and — 2J, respec- 
tively. For each induction direction, the signs of the Jb 
and Jd contributions to b\ c ^ + 2J and B X ^p + J are the 
same, whereas the signs of the J a and J e contributions 
to B l 2^ + 2J and B^p + J are the opposite. 

In Figs. 6-10, we plot M/j versus jB/\J\ for five low- 
T cases of AFM s x = 1 dimers, taking k B T/\J\ = 0.03. 
These curves all exhibit the universal features predicted 
by Eqs. jlB} and ijlijl. The corresponding Cv/kB versus 
-fB/\J\ curves, which also obey Eqs. (|4"0"|l - (f4"2")l . are pre- 
sented elsewhere, [ijj In these figures, only one of the five 
anisotropy interactions Jj is non- vanishing, and we take 
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FIG. 6: Plot at k B T/\J\ = 0.03 and J b /J = 0.1 of M/7 
versus r yB/\J\ for the AFM spin 1 dimer, with the same curve 
notation as in Fig. 2. 

Jj/J = 0.1, for j = b,d,a,e, and c, respectively. Unlike 
the case of s± = 1/2 dimers, for which J c merely renor- 
malizes J<j, for s\ = 1 each of these interactions leads to 
distinct anisotropy effects in the low-T thermodynamics 
arising from the B\ c 1 (9, <f>). 

We first show the results for the global anisotropy in- 
teractions. In Fig. 6, we plot M{B) for Jt, = 0.1J (and 
all other Jj = 0) at k B T/\J\ = 0.03 for the AFM s x = 1 
dimer. B\ c 1 (9, 0) is largest and increases monotonically 
with level-crossing number s for 9 — 0, but is indepen- 
dent of s for 9 = 7r/2, nearly quantitatively consistent 
with Eqs. (|59|) - (|62[1 . In Fig. 7, the corresponding re- 
sults for Jd = 0.1 J (and the other Jj = 0) are shown. 
Again, B^^O, 0) is largest and increases monotonically 
for 9 = 7r/2, and is nearly independent of s for 9 = 0, 
as for s = 1/2, nearly quantitatively consistent with Eqs. 
dnnj)-(in2J • Note that -B^i(0, 0) is nearly indistinguishable 
from the isotropic case. 

Next, in Figs. 8 and 9, we show the corresponding 
curves for M (B) at low T for the local anisotropics J a — 
0.1J and J e = 0.1J, respectively, with the other Jj — 
0. In contrast to the global anisotropies, B^-^, 0) + sJ 
changes sign for both 9 — 0,7r/2, nearly quantitatively 
consistent with Eqs. $5$$-^ ■ We note that 5^(0, 0) is 
nearly indistinguishable from the isotropic case for J e =/= 
shown in Fig. 9. 

Finally, in Fig. 10 , we show the corresponding M(B) 
curves for J c = 0.1J. With this combination, B\ c 1 (9, 0) 
is independent of s for both 9 = 0, 7r/2, but is larger 
for 9 = 7r/2, nearly quantitatively consistent with Eqs. 
(|59ll - ljr]2*)l with J d , J e — * ± J c /2, respectively. 

In short, the case Si = 1 is sufficient to exhibit the 
very different behaviors obtained from the single-ion, lo- 
cal spin anisotropy interactions from those obtained from 




9H B B/|J| 

FIG. 7: Plot at k B T/\J\ = 0.03 and J d /J = 0.1 of M/7 
versus 'yB/\J\ for the AFM spin 1 dimer, with the same curve 
notation as in Fig. 3. 




gn B B/|J| 

FIG. 8: Plot at k B T/\J\ = 0.03 and J a /J = 0.1 of M/7 versus 
7-B/l J\ for the AFM spin 1 dimer. The curve notation is the 
same as in Fig. 2, except that the isotropic case (dotted) 
curve is for J a — 0. 

the global spin anisotropy interactions. As si increases 
beyond 1, the situation becomes not only more compli- 
cated, but also more interesting, as shown in the follow- 
ing. 

VI. EXACT NUMERICAL RESULTS FOR SPIN 

5/2 

For si = 5/2, one of the cases of greatest experimental 
interest, when Ti. a and Ti. e are present, none of the allowed 
s,m values is a true quantum number. That is, Ti. a and 
H e cause all of the states with nominally odd or even s to 
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J e /J=0.1, k Q T/|J|=0.03 
AFM S=1 

;B||x 



1.6 2.0 

g^ B B/|J| 



— i 

2.4 



J b /J=0.1, k B T/J=0.03, S=5/2dimer 
AFM 




2 3 4 

gjx B B/|J| 



FIG. 9: Plot at k B T/\J\ = 0.03 and J e /J = 0.1 of M/7 versus 
7-B/l J| for the AFM spin 1 dimer. The curve notation is the 
same as in Fig. 2, except that the isotropic case (dotted) 
curve is for J e = 0. 



HQ 



J/J=0.1, k B T/|J|=0.03 



AFM S=1 





0.8 



1.2 1.6 2.0 

g^B/|J| 



FIG. 10: Plot at k B T/\J\ = 0.03 and J c /J = 0.1 of M/7 
versus 7-B/l J| for the AFM spin 1 dimer. The curve notation 
is the same as in Fig. 2, except that the isotropic case (dotted) 
curve is for J c — 0. 



mix with one another. For B\\i for i = x,y, z, this sim- 
plifies in the crystal representation, as for s\ = 1, since 
only states with odd or even m values in the appropri- 
ately chosen representation can mix. By using symbolic 
manipulation software, it is possible to solve for the exact 
eigenvalues of the s\ = 5/2 dimer. However, because the 
analytic expressions for the eigenvalues are much more 
complicated than those for s± — 1 presented in Appendix 
A, we shall not attempt to present them, but will instead 
focus upon their numerical evaluation for specific cases. 
To first order in the Jj, the first three level crossings 



FIG. 11: Plot at k B T/\J\ = 0.03 of M/7 versus -)B/\J\ for 
the AFM spin 5/2 dimer with Jb/J = 0.1. The curve notation 
is the same as in Fig. 2. 
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FIG. 12: Plot of M/7 versus jB/\J\ for the AFM spin 5/2 
dimer at k B T/\J\ — 0.03 with J d j J — 0.1. The curve notation 
is the same as in Fig. 3. 



for B\\i with i — x,y,z are 
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FIG. 13: Plot at k B T/\J\ = 0.03 and J a /J = 0.1 of M/7 ver- 
sus 7-B/| J| for the AFM spin 5/2 dimer. The curve notation 
is the same as in Fig. 8. 



FIG. 14: Plot at k B T/\J\ = 0.03 and J c /J = 0.1 of M/7 ver- 
sus 7-B/| J| for the AFM spin 5/2 dimer. The curve notation 
is the same as in Fig. 10. 



In Figs. 11-14, we plot M/7 and C V (B) versus jB/\ J\ 
for four low-T cases of AFM s\ = 5/2 dimers, J& = 0.1 J, 
Jd = 0.1J, J a = 0.1J, and J c = 0.1J, respectively, and 
the other Jj = 0, taking k B T/\J\ = 0.03. Each of these 
curves exhibit the universal step behavior predicted in 
Eqs. (|4"5|) and fPtJl . Corresponding Cy/fce versus r )B/\J\ 
curves also exhibit the universal double peak behaviors 
predicted in Eqs. (|40() - (|42|l . and are shown elsewhere. [2^ 
In Fig. 15, examples of M(9) at fixed B and = 
are shown. Figures 13-15 are sufficient to distinguish 
the more interesting local spin anisotropy effects in AFM 
dimers with higher s\ values from the non-existent or less 
interesting ones present with S\ = 1/2,1, respectively. 
In Figs. 11-14, the solid and dashed curves represent the 
B\\z and B\\x cases, and the dotted curve represents the 
isotropic case, Jj = OVj, as in Figs. 2-10. 

We first examine the global spin anisotropy effects of 
Tit, and Tid in Figs. 11 and 12. These figures exhibit the 
same behavior shown for s\ = 1/2, 1 in the corresponding 
Figs. 2, 6 and 3, 7. Note that B\ c 5/2 {9,0) is largest for 
9 = tt/2 with J b = 0.1J and for 6 = with J d = 0.1J, 
and increases monotonically with s, as for s% = 1 in both 
cases, nearly quantitatively consistent with Eqs. 1631) - 
(P|l. By contrast, B^ 52 (6, 0) for 6 = tt/2 with J b = 0.1J 
and for 9 = for Jd = 0.1J are nearly indistinguish- 
able from the isotropic case, also nearly quantitatively 
consistent with Eqs. (|fj3^l - lltjS|l . 

In contrast, the local field anisotropy interactions show 
very different and much more interesting behaviors. In 
Fig. 13, we present our results for the effects of the 
single- ion axial anisotropy interaction 7i Q , Eq. (J5J. As 
in Fig. 8 for si = 1, B^ 5 , 2 (6, 0) 4- sJ changes sign 
with increasing s. B l f 5/2 (0,0) + J < 0, whereas for 



s > 2, _B' C 5 ^ 2 (0, 0) + sJ > and increases monotoni- 
cally with s. For 9 = tt/2, nearly the opposite situa- 
tion occurs. i?' c 5 y 2 (7r/2, 0) + sJ is positive for s = 1 
and decreases monotonically with increasing s. In both 
cases,|B^ c 5 ^ 2 (0, 0) -I- sJ\ is a minimum for s = 2. These 
local axial anisotropy effects, consistent with Eqs. (1031) - 
(|68fl . are very different than the global anisotropy ones 
pictured in Figs. 11 and 12. They are also much richer 
and interesting than the corresponding case for si = 1 
pictured in Fig. 8. 

In Fig. 14, we present our M/7 versus ^B/\J\ re- 
sults for the case of the local azimuthally anisotropic 
exchange interaction, H c , Eq. ©, evaluated for AFM 
dimers at ksT/\J\ = 0.03 with J c = 0.1J and the re- 
maining Jj = 0. B^l for this case is obtained from 
Eqs. JHSJ-lEHl by setting Jd,J e = ±J c /2, respectively, 
and J a = Jb = 0. B l ^ 5 / 2 (0, 0) is nearly indistinguish- 
able from the isotropic case, as in Fig. 10 for si = 1, 
except for some minor curve shape effects far from the 
step midpoints. B^ 5 ^ 2 (tt/2, 0) + sJ is always positive, 
as for si = 1 shown in Fig. 10, but has a minimum at 
s = 3. This is also in stark contrast to the monotonic 
global anisotropy behavior for si = 5/2 seen in Figs. 11 
and 12. Both of these behaviors are nearly quantitatively 
consistent with Eqs. (|63fl - H68[) as modified to include J c . 

Although not pictured explicitly, the behavior for J e = 
0.1J with the other Jj = is rather like that of the J a = 
0.1J curves pictured in Fig. 13 with z <-> x, differing in 
ways similar to those differences between the J a = 0.1 J 
and J e = 0.1J curves pictured for s% = 1 in Figs. 8 and 9. 
As indicated in Eqs. (ESJ-lEHl), tne ^5/2^^ °) are nearl y 
independent of J e . However, the B l ^ 5 ^ 2 (n/2, 0) + sJ with 
J e = 0.1J are nearly three times as large as for J a = 0.1 J 
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B. 

In these rotated coordinates, the Zeeman interaction 
—bS z >, where b = -fB, is diagonal. We therefore denote 
this representation as the induction representation. The 
Hamiltonian Ti' in this representation is given in Ap- 
pendix B. In the induction representation, we choose the 
quantum states to be |<£™)- In the absence of the four 
anisotropy interactions Jj, Ti.' — Ti' is diagonal, 



Ha 



where 



£™' (0) = -Js(s + l)/2 -mb. 



(69) 



(70) 



FIG. 15: Plot of M/7 at 4> = versus 9/iv near the sec- 
ond step at jB/\J\ = 2.0,2.1 and k B T/\J\ = 0.03 for each 
Jj/J = 0.1 with i 
j, 7B/I J| values. 



case, so that the minimum \B l /r 5 / 2 (ir/2, 0) + sJ\ is also a 
minimum for s = 2. More details are given in Section 
VII and Appendix D. 

In addition, the angular dependencies of M are differ- 
ent for each of the four Jj we presented for si = 5/2. 
In Fig. 15, we present the results for M(B,9,(f) — Q)/j 
versus 9/ir near the second level crossing at r yB/\J\ — 
2.0, 2.1 and k B T/\J\ = 0.03 for each of the four si = 5/2 
AFM magnetization cases pictured in Figs. 11-14. Note 
that M(B, tt - 0, 0) = M (B, 9, 0). The J a and J b curves, 
while rather similar at jB/\J\ = 2.0, are very different 
at jB/\ J\ = 2.1. Hence, M(B) depends strongly upon 
the particular type of spin anisotropy. 



VII. ANALYTIC RESULTS FOR WEAKLY 
ANISOTROPIC DIMERS OF ARBITRARY SPIN 

A. Induction representation eigenstates first order in 
the anisotropics 

Since the diagonalization of the Hamiltonian matrix 
is difficult for an arbitrary magnetic field direction com- 
bined with an arbitrary combination of spin anisotropy 
interactions, and must be done separately for each value 
of si, it is useful to consider a perturbation in the rel- 
ative strengths Jj/J of the anisotropy interactions. We 
nominally assume \ Jj/J\ <C 1 for j — a,b,d,e. However, 
to compare with low-T magnetization and specific heat 
experiments at various applied field directions and mag- 
nitudes, one cannot take the magnitude of the magnetic 
induction to be small. In order to incorporate the mag- 
netic induction of an arbitrary strength and direction ac- 
curately, we therefore rotate the crystal axes (x, y, z) to 
(x' , y', £'), so that B = Bz' . The rotation matrix and a 
brief discussion of its ramifications are given in Appendix 



The operations of the remaining terms in Ti.' on the 
eigenstates are given in Appendix C. The first order 
correction to the energy in the induction representation, 
£ s m <« = (<p?\K'\<pY), is found to be 

= -y [2a(. + l)-l]-£[.(- + l)-l] 

Js,si 

+ J ^-[m 2 + s(s + l)-l] 
+ I[ s ( s + l)-3m 2 ] 

X (Jb',a COs2 8 + Jdle 1 ^ 9 c °s( 2 ^)) , (71) 

where 



J b,a 

Js,S! 

J d,e 



= Jb 
= Jd 



si Ja j 



(72) 
(73) 



arises from the 



and ot s S1 is given by Eq. I|29|) . 

Since the 9, dependence of E™'^ 1 ' 
term proportional to J^ 1 cos 2 9 + sin z 9 cos(20), 
it is tempting to think that the thermodynamics with 
Jde 1 = and B\\z are equivalent to those with J^ 1 = 
and _B||i;. However, as shown explicitly in the following, 
the 9, (^-independent parts of Eq. l(7T|) strongly break this 
apparent equivalence, causing the Bl c Si (9, (j>) for those 
two cases to differ. This implies that Jb and Jd are in- 
equivalent, as are J a and J e , even to first order in the 
anisotropy strengths. 



B. First order thermodynamics 

In the Hartree approximation, the anisotropy interac- 
tions are included to first order only. In this approxima- 
tion, s and m are still good quantum numbers, so the 
partition function 



2,S! 



zil) = E E 



-0E" 



(74) 



s—0 m— - 
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where E™ Sl = Ef m + E™s{ 1] '. Although it is difficult to 
perform the summation over the m values analytically, it 
is nevertheless elementary to evaluate Z numerically for 
an arbitrary B, 8, <fi, and T from the eigenstate energies. 
The magnetization in the Hartree approximation is 

M«(B,M) = ^ryE E rne-^. (75) 

s— m=—s 

Similarly, the specific heat in the Hartree approxima- 
tion is 

C$\B,9,*) « T^^E E (EZJ 2 e-^ 
\^ { ') s =0m=-s 

-(E E C.e^"') 2 ]- (76) 

s— m— — s 

As a test of the accuracy of this Hartree calculation, we 
have compared the Hartree and exact M(B) obtained for 
the si = 5/2 dimer with = 0.1J and B\\z at various 
T values in Fig. 1161 The corresponding comparison be- 
tween the Hartree and exact Cy{B) is shown in Fig. El 
We see that the curves evaluated using the Hartree and 
the exact expressions for M and Cy with s\ — 5/2 are 
indistinguishable at ksT/\J\ = 0.03. The GV curves are 
noticeably different at k B T/\J\ = 0.1 for >yB/\ J\ < 0.4, 
and at ksT/\J\ — 0.3 they are noticeably different for 
jB/\J\ < 2.6. Corresponding noticeable differences in 
the M curves at the same B values appear at T values 
roughly three times as high as in the Cy curves. 

At very low T, ksT/\J\ <C 1, the most important 
states in this perturbative scheme are the minima for 
each s value, Eg a , which determine the level crossings 
in the Hartree approximation. As T — * 0, we can ig- 
nore all of the m ^ s states in Eqs. JSJ-JZEJ- This 
two-level approximation is the basis for the universal be- 
havior given by Eqs. (|40|l - i|44|) . which fits all of the exact 
curves we presented reasonably well. 




gn B B/|J| 

FIG. 16: Comparison of M/y versus jB/\J\ obtained using 
the Hartree asymptotic form (dotted) with the exact calcu- 
lation (solid), for the si = 5/2 AFM dimer with J d = 0.1J, 
J a = Jb = Je = 0, at k B T/\J\ = 0.03,0.1,0.3,0.8, as indi- 
cated. Inset: expanded view of the region < jB/\J\ < 1.5. 
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C. First order response functions 

Inelastic neutron scattering cross-section 



FIG. 17: Comparison of Cv/kB versus jB/\J\ obtained using 
the Hartree asymptotic form (dotted) with the exact calcu- 
lation (solid), for the si = 5/2 AFM dimer with J d = 0.1J, 
Ja = Jb = Je = 0, at k B T/\J\ = 0.03,0.1,0.3, as indicated. 



Two response functions relevant for the study of 
SMM's are the inelastic neutron scattering cross-section 
S(B,q,u>) and the electron paramagnetic resonance 
(EPR) susceptibility \{B,lo) in strong magnetic in- 
ductions B. The low-T inelastic neutron cross-section 
S(B, q,u>) can be evaluated from 

3 

S(B,q,u) = E OW-Ma) / 

a, 0=1 71 

xTr(e-P u S a (q,t)S (q,O)), (77) 
S a (q,0) = S 1 ^ d + S 2 ^ a e^ d , (78) 



where 2d is the vector between the spins of a dimer and q a 
is the ath component of q/q. It is customary to write the 
operators in the crystal representation, for which d = dz. 
[lj,|2y| With inelastic neutron scattering, one can probe 
the dimer with a strong magnetic field at various direc- 
tions with respect to both the dimer axis and the scat- 
tering plane. To the extent that the eigenstates {|0 n )} 
and their energies e„ can be evaluated exactly, 

3 ™-i 

S(B,q,u) = E (<W - fe<Z/3) E e_/3e " 

ex,/3=l n,n' — 1 
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x6(u + e„ - e„/) 

x (<f> n | S a (q, 0) \<j> n ,) (<t> n , | Sj(g,0) |<M, 

(79) 

5 a (g,0) = t75 a (q,0)^t, (80) 

where U is the unitary operator that diagonalized TL. 
When exact expressions for S(B, q, ui) are tedious to ob- 
tain, it is straightforward to obtain it in the Hartree 
approximation, S^(B,q,uj). But to do so, it is eas- 
iest to define the axes in the induction representation, 
for which {|0 n )} = {|<^™)} an( i the e Q = x',y',z' for 
q = 1,2,3, respectively. In the Hartree approximation 
obtained by setting the |y>™) = \ip™), the bare wave func- 
tions, (B, q, ui) for an arbitrary B, 9, <f> is then 

2si s 
s— m=—s 

x (cos 2 (g • d) sin 2 6>6, g ^^ (0) (w, 0, 0) 

+ sin 2 (q • d) sin 2 fl&.g.F™^ 05 (w, 0, 0) 

2/ ,n2 — sin 2 b9 m , s ( ), . 
-cos z (<j-d) ^^ 3 , Sl (u,0,(l>) 



+ sin 2 (q • d) 



4 

2 - sin 2 i 



b,q ^-m,s(0) 



4,si 



(81) 



where q ■ d = qd cos # g is invariant under the rotation, 2d 
is the vector separating the dimer spins, 9b lQ is the angle 
between q and B, and 



= m 2 S(to), (82) 

^; (0) = E 
<j'=±i 

^ (0) = e 5(^+^ m sl -^r)K m ) 2 ) (84) 

er = ±l 

^; (0) - E *(«+*Ci 



<r,o-'=±l 



- (^:++%tr' m j . (86) 

where the coefficients are given respectively by Eqs. (I21II . 
JTJ, and H20fl . and we have suppressed the uj,8,(j) argu- 
ments of the J-'^sl^ functions for simplicity of presen- 
tation. In Eq. (|81|l . the only dependencies upon the 
anisotropy energies is in the first-order eigenstate ener- 
gies E™ S1 = ET m +ET^ ] given by Eqs. EJ9 and (EJ), 
respectively. 

As described in detail in Appendix C, corrections to 
the wave functions first order in the anisotropy inter- 
actions lead to the "extended Hartree" approximation, 
S^ lc \B, q, w), which exhibits additional transitions with 



strengths first and second order in the anisotropy ener- 
gies. S( lc ) (B, q, to) contains terms arising from the trace- 
less diagonal and off-diagonal parts of the a, (3 component 
matrix, and from the contributions of the diagonal part 
of the component matrix containing corrections to its 
trace. These latter terms are proportional to the same 
four functions of q as in Eq. (|81|l . and can be sepa- 
rated experimentally from the other terms by choosing 
the scattering plane such that q _L B (or sin0j, )9 = 1) 
for all {9, <f>) and by averaging over all azimuthal angles 
4>b.q- Then, the transitions from (s,m) to (s',m'), where 
s' = s ± 2, 3 and m' = m, m ± 1, m ± 2, m ± 3 with lead- 
ing amplitudes second order in the anisotropy interac- 
tions upon depend upon the the local anisotropy interac- 
tions J a and J e through the three functions / n (0, 0) with 
n = 3, 6, 7 given in Appendix D. Hence, if these transi- 
tions can be identified, a study of the (9, (f>) dependen- 
cies of their strengths can provide direct measurements 
of those parameters. Then, studying the weak transi- 
tions (s,m) (s',m') with s' = s, m± 1,±2, ±3 and 
s' = s±l,m' = m±2 can provide further measurements 
of the global anisotropy interactions Jd and Jd- 



2. Electron paramagnetic resonance susceptibility 

In an EPR experiment, one applies a strong mag- 
netic field, plus a weak oscillatory transverse field, lead- 
ing to the overall induction B = Bz! + B±[x' cos(ut) — 
ay' sin(wi)], where we assume the oscillatory induction 
precesses clockwise (counterclockwise) for a = ±1.[26) 
For a weak transverse induction Bj_, one measures the 
resulting linear response X—o-,a(B , u>) given by 



?7 

JUt 



dt 
2V 



■ i( - a ' +i,, ' )t Tr(e- f} ' H [S- a (t),S a (0)} 



-iUt 



(86) 
(87) 

where ?y = 0+, the global spin raising and lowering op- 
erators S± are defined in the induction representation, 
as described in Appendix B and Z is the partition func- 
tion, and [a, b] = ab— ba is the commutator. In principle, 
an exact expression for X—a-,a(B,u>) can be obtained by 
starting in the induction representation with the basis 
diagonalizing 7i with the unitary operator V, 
VHV^ = H', and the new basis {|0n)} is obtained from 
\4>n) = leading to 7i'|</>„) = e n \<j> n )- One then has 



X- 



(B,lj) = -h= 



7 



y 

irZ V 

n , n' — 1 



X (0„ | (o) | <M (<M s CT (o) 1 0* 
i 



- 17] 



(88) 



13 



where S a (0) = VS <J (0)V" <! . It is then elementary to ob- 
tain X-a ai 10 ) m the Hartree approximation valid to first 
order in the anisotropy parameters, which is well-behaved 
for all induction values. The imaginary part of X—a,o is 
then 

s— m— — s 

[s(s + 1) — m 2 — am] 

-5(E™Y-E™ Si +oj)), (89) 

where E™ s 1 = ET m + E'T^ ] is given by Eqs. (25 and 
(|71|l and Z^> is given by Eq. 1)74)) . The arguments of the 
(5-function give rise to the resonant frequencies first order 
in the anisotropy energies, or to the first order resonant 
magnetic inductions, 

lB w = ±UJ+ (^±^lg ssi{eA) , (go) 

g., ai (0,4>) = J b % Sl (l-3cos 2 0) 

-3Jjf sin 2 6»cos(20), (91) 

where J^ 1 and J^ 1 are given by Eqs. ({?2*jl and 
respectively, and we have ignored the processes that give 
rise to finite transition widths. [27l| By varying the direc- 
tion and magnitude of B, it is possible to obtain sufficient 
information to fit all of the parameters in the model. At 
low T the m — s states dominate, so that only a = — 1 
is allowed. By increasing B past the sth level crossing, 
one can probe the effective sth state of the AFM dimer 
using EPR. 

In Appendix C, we extend these results to include 
the corrections to the wave functions first order in the 
anisotropy interactions. With these corrections, addi- 
tional EPR transitions are present, which can provide 
additional information useful in experimental identifica- 
tion of the anisotropy interactions. It is shown that there 
are 10 additional resonant magnetic induction strengths 
of the forms 

jB^l = ±a n u> + b n g s , Sl {9,4>) 

+c n h s , Sl (9,<t>), (92) 

= -(J+2J b + J a ) 

+J b % Sl (l + cos 2 #) 

+J s d f sin 2 6>cos(2</>), (93) 

where the a n ,b n , and c n are listed in Appendix C. We 
note that these additional resonant inductions are first 
order in the anisotropy interactions, but have amplitudes 
that are second order in the anisotropy interactions, so 
that in most cases, they may be difficult to detect, but 
their detection in certain materials would provide a clear 
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FIG. 18: Plot of the exact eigenvalues Cn/\J\ versus jB/\J\ 
for st = 1 AFM dimers with J a /J = 0.1 and 6 = tt/4. The 
arrows represent the strongest EPR transitions at their cor- 
responding field strengths. 

signal of the presence of significantly strong anisotropy 
interactions. In addition, resonant inductions with c n ^ 
are generally very large. 

To illustrate one example of an effect of local spin 
anisotropy upon the EPR transitions, we consider the 
simple case of an si = 1 AFM dimer with the only non- 
vanishing anisotropy energy J a /J = 0.1 and 9 = tt/4. 
The energy levels for this system are pictured in Fig. 18. 
In Fig. 19, the EPR transition energies versus "fB/\J\ 
for the sx = 1 AFM dimer with J a /J = 0.1 and 9 = tt/A 
are shown. The widths of the lines are proportional to 
the strengths of the EPR matrix elements. Note that the 
ground state at B = is approximately \ip$), but its en- 
ergy is slightly negative due to an admixture of this state 
with the nominal IV 1 ™) states with even m, as shown in 
Appendix A. The EPR matrix elements are very small 
at these B values, as the transitions only exist due to the 
admixture of the wave functions. At higher B strengths, 
the ground state goes through two level crossings, with 
the first level crossing being to the lowest energy s = 1 
state, which is nominally \4>l), with some mixture of the 
other {ip™) states. The leading transition is to the nomi- 
nal \ipi) state. Then, the second level crossing causes the 
ground state to be the nominal |-0 2 ) state, which is of 
course modified by the mixing with the nominal |^>q) and 
the other nominal |?/>™) for m = 0, —2. Note that near to 
jB/\J\ = 3, there is a level repulsion. The leading EPR 
transition is to the nominal \ip\) state, and because of 
the strong matrix elements, the splitting of the energies 
due to the level repulsion should be observable in EPR 
experiments. 

D. Level crossings first order in the anisotropy 
energies 

We can find an expression for the s th AFM level cross- 
ing at the induction B^ii to first order in the anisotropy 
interactions for a general s\ spin dimer by equating 
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FIG. 19: Plot of the EPR transition energies (e n — e\)/\J\ 
versus ~]B/\J\ from the ground state for the si = 1 AFM 
dimer with J a /J = 0.1 and 6 = 7r/4. The widths of the lines 
are proportional to the strengths of the matrix elements for 
the transitions. 



E 



s,(0) 



E 



to E. 



-i,(o) 



-i,(i) • u- 

-i,a\ . yielding 



(4s - 3) 



[Jb cos 2 9 + Jd sin 2 9 cos(20)] 



+3c s , Sl [J a cos 2 9 + J e sin 2 9 cos(20)], (94) 
[3 + 3s - 5s 2 - 4s 3 + 4s x (si + 1)] 



2(2s+ l)(2s + 3) 



(95) 



This expression is consistent with those obtained for 
si = 1/2,1,5/2 given by Eqs. JSTJl, and (gg)- 

ffBfy. in Figs. 11-14. We note that -yS^i^ contains the 
0, (/(-independent terms, — Js — Jb/2 — c s , Sl J a , which dis- 
tinguish Jb from Jd and J a from J e . 

In particular, we note that the single-ion anisotropy 
interactions behave very differently with increasing step 
number than do the global anisotropy interactions, espe- 
cially for large Si. For the three cases we studied in de- 
tail, for si = 1/2, ci .1/2 = 0, so that the local anisotropy 
terms are irrelevant, for si = 1, ci t i = | and C2.1 = — | 
have different signs, and for si = 5/2 as in Fe2 dimers, the 
first three c Si 5/2 coefficients are y|, — and — 1|, respec- 
tively, the second being an order of magnitude smaller 
than the other two, and opposite in sign from the first. 
For si = 9/2 dimers such as [MJI4J2, the first four c Si g/ 2 
are |, — |, and — 1|, which changes sign between 
s = 2 and s = 3, where its magnitude is a minimum. 

This is in sharp contrast to the global anisotropy in- 
teractions, for which the analogous coefficient (4s — 3)/2 
increases monotonically with s, independent of si. These 
differences should be possible to verify experimentally in 
careful low-T experiments at high magnetic fields applied 
at various directions on single crystals of those si = 5/2 
Fe2 and si = 9/2 [Mn,^ dimers for which |J| is suffi- 
ciently small. 



E. Level crossings to second order in the anisotropy 
energies 



To aid in the analysis of experimental data, we have 
extended this perturbative calculation to second order 
in each of the four principle Jj. Since we expect the 
\Jj/J\ < 0-1 in most circumstances, this extension ought 
to be sufficient to accurately analyze most experimen- 
tally important samples. Since the c SjSl magnitudes can 
become larger than unity for si — 5/2 and especially for 
81 = 9/2, we expect second order corrections to be sig- 
nificant in those cases. To second order in the anisotropy 
interactions, the eigenstate energies _E™sP^ are given in 
Appendix C. We note that the E^sf 1 contain divergences 
at jB/\J\ = 0,2s- 1,2s + 3, s- 1/2, and s + 3/2, so that 
near to those values, one would need to modify the per- 
turbation expansion to take proper account of the degen- 
eracies. Hence, the expressions for E™s? cannot be used 
in the asymptotic expressions for the thermodynamics, 
Eqs. (HlJ-(|2Il!- However, as the sth AFM level crossing 
occurs approximately at ~{B/\J\ = s, which is far from 
any divergences, we can safely use this second order ex- 
pansion to obtain an expression for the level crossings 
second order in the anisotropy interaction energies. We 
find 



IB 1 *? 



™,(2) _ F s-U2) 



B=-Js/ 1 



(96) 



where 7B^' is given by Eq. (|94J) . 

The full expression for the second term in Eq. (|96|l is 
given in Appendix D. From this expression, it is easy to 



see that for s x = 1/2, £^ 2) 2 = ^f[ 2) ( 



(2) (2) 

where f{ and are given in Appendix C. For B\\z, 



E 



lc(2) 
1,1/2 



iJ 2 J\J\, and with B\\x,y, E\ c [% = -\(J b ± 



/2 



Jd) /\J\i i n agreement with the expansion to second or- 
der in the Jj of our exact formulas in Eq. I)57|l . However, 
the second order functions have more complicated 9, </> 
dependencies than do the first order in Eq. 

We have also explicitly checked each formula in Appendix 
D for si = 1 and s = 1,2. Thus, Eq. I® is a highly 
accurate expression for the full 9, 4> dependencies of all 
s = 1, . . . , 2si + 1 level crossings of a single crystal dimer 
of single ion spin s x . 

By superposing this non-universal level-crossing for- 
mula, Eq. (|96(l . combined with the universal behavior 
presented in Eqs. I|4Q(I - I|44(I . it is easy to use our results 
in accurate fits to experimental data at low temperatures 
and high magnetic induction strengths. 
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VIII. AFM DIMERS WITH STRONG 
ANISOTROPY INTERACTIONS 

Finally, we consider some cases of strong anisotropy 
interactions, in which one or more of the Jj is compa- 
rable to J in magnitude. The cases of interest are those 
for which the anisotropy interactions can remove the level 
crossing effects in AFM dimers that give rise to the 2si + l 
M(B) steps and Cv(B) double peaks. These cases gen- 
erally occur for strong FM anisotropy interactions, but 
there are some examples in which they can occur with 
strong AFM anisotropy interactions. We first consider 
the case of s\ = 1/2 dimers, for which the situation can 
be analyzed analytically, at least for B\\i. Then, we con- 
sider the si = 5/2 case numerically. 




9H R B/|J| 



A. Analytic and numerical results for si = 1/2 dimers 

For the s\ = l/2 dimer, the energies for B||i are given 
by Eqs. i)45 [l -l]51 [l . We are interested in examining the 
cases in which the anisotropy is strong enough to cause 
the level crossing to disappear. Strong anisotropy with 
Jb and Jd having the same sign as J do not differ sig- 
nificantly from the weak coupling cases, as they do not 
remove the level crossing, and do not affect significantly 
the heights of the magnetization step and the number 
and shapes of the C\r{B) double peaks, but just shift 
their positions, as for weak anisotropy. There are then 
three AFM (J < 0) cases of interest. These are: (1) 
Jb = Jd = \J\/2, (2) J b = |J|, and (3) J d = \J\. The 
M(B) and Cy(B) curves for these cases with B\\z at 
the very low-T value ksT/\J\ = 0.01 are shown as the 
solid, dotted, and dashed curves in Fig. 20 and 21, re- 
spectively. 

For case (1), J, = Jd = \J\/2, the same four eigenval- 
ues are obtained for B\\z and Assuming J < 0, 
these are ranked from highest to lowest as 



£3 
ei 



-J/2 + v/6 2 + J 2 /4, 
J, 



0. 



-J/2 - v/6 2 + J 2 /4. 



(97) 
(98) 
(99) 
(100) 



Note that t\ — €2 only at b = (where 63 = 64), t\ de- 
creases with increasing b, and 64 increases with increasing 
b, so the levels get further apart with increasing b. In this 
case, M(B) shown as the dotted curve in Fig. 20 is broad, 
even at the very low T value plotted. In addition, Cy (B), 
shown as the dotted curve in Fig. 19, exhibits a single 
Schottky-like anomaly, but with a different shape than 
that of the standard Schottky case, which arises from a 
splitting of the two lowest energies linear in b at b = 0. In 
this case, the peak value of Cy/ks is 0.439229, the same 
as the uniform double peak value, but the peak position 
is as b = x /2c|J|//3, where c = 1.19967864. 



FIG. 20: Plots of M/7 versus jB/\J\ at k B T/\J\ = 0.01 with 
J b = \ J\ (solid) for both B\\x,z, J b = J d = |J|/2 (dotted), 
and J d = I J\ with B || i. (dashed) for the si = 1/2 AFM dimer. 
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FIG. 21: Plots of Cv/k B versus jB/\J\ at k B T/\J\ = 0.01 
with J b = \J\ (solid) for both B\\x, z, J b = J d = | J|/2 (dot- 
ted), and J d = \ J\ with B\\z (dashed) for the si = 1/2 AFM 
dimer. 



Next, we examine case (2), J, = |J|, Jd = 0. In this 
case the energies are different for the two field directions. 
For Bllz that 



Q 
£3 
£2 



J. 



b. 

0. 



(101) 
(102) 
(103) 
(104) 



At b = 0, these energies comprise a triply degenerate 
(or spin 1) ground state at and an excited +| J| state. 
With increasing 6, t\ decreases linearly, and £4 increases 
linearly, crossing 64 at b = | J|, which does not affect the 
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low-T thermodynamics. In this case, M(B) is a sharp 
step at B = 0, shown as a solid curve in Fig. 20, and 
Cy(B) is a conventional spin 1 Schottky anomaly at low 
b, as shown by the left solid curve in Fig. 21. From Eqs. 
(|TU5|l - (jTU3|) . C v /k b has a maximum value at 0.637203 at 
b = 1.8806775//3. Note that the larger maximum value 
than in the other Cy(B) cases arises from the higher 
degeneracy as b — > 0. 

For b||x, the energies for case (2) are 



£4 
£2 



-J/2 + v^ 2 + J*/*, 



= 0, 



= -J/2 - v/6 2 + J 2 /4. 



(105) 
(106) 
(107) 
(108) 



At 6 = 0, the ground state is again triply degenerate with 
energy 0, and the excited state has energy | J|. With in- 
creasing b, t\ separates from £2 and £3, curving below 
them, and £4 increases quadratically. For this induction 
direction there is no level crossing, but there is some b 
anisotropy at finite T, because of the quadratic versus 
linear b dependencies. Cy/ks has the maximum value 
0.761802 at b = ^2 Cl \J\/(J, where c x = 2.654658, which 
is noticeably different from the case with B\\z. The 
M(B) and Cy{B) curves are the right solid curves in 
Figs. 20 and 21, respectively. 

Now we consider case (3), Jd = — J, Jb = 0. Again, 
there are slight differences for B||z and B||x. For B||z, 
the eigenstate energies are 



£3 



-J+y/P + J 2 , 

-J 

0, 



-J- Vb 2 + J 2 . 



(109) 
(110) 
(111) 
(112) 



At b — 0, the ground state is doubly degenerate at en- 
ergy 0. The other two states have energies | J| and 2| J|, 
respectively. As b increases, t\ decreases below £2, and £4 
increases monotonically. Hence, there is no level cross- 
ing, but at low T, two levels are relevant, in a fashion 
slightly different from a standard Schottky specific heat 
anomaly. The low-T M(B) and Cy(B) curves for this in- 
duction direction are shown as the dashed curves in Figs. 
20 and 21, respectively. We note that they are rather sim- 
ilar from those of case (1), but have the broadest M(B) 
increase and Cy (B) peak of the three cases considered. 
Cv/ks has the peak value 0.439229 which is the same 
as the uniform double-peak value, but its position is at 
b = 2y/c\J\//3, where c = 0.19967864. 
For Bllx, we have 



£ 4 = -2J, 

£3 = -J/2 + Vb 2 + J 2 /4, 

(-2 = 0, 

£1 = -J/2 - ^b 2 + J 2 /4. 



(113) 
(114) 
(115) 
(116) 
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FIG. 22: Plots of M/7 versus -yB/\J\ at k B T/\J\ = 0.01 with 
Jb = Jd = c\.J\ for the AFM si = 5/2 dimer. The cases 
c = 0.1 (solid), c = 0.2(dashed), c = 0.3 (dotted), c = 0.4 
(dot-dashed), c = 0.45 (short dashed), and c = 0.5 (short 
dotted) are shown. 



At 6 = 0, the eigenstates are the same as for the field in 
the z direction, with a doubly degenerate ground state at 
0, and excited states at | J| and 2| J|. However, the field 
dependence is a bit different than for the z direction. 
Again the ground state decreases with increasing b, but 
the energy scale of the curvature is | J | / 2 instead of | J | . In 
addition, £3 can cross £4 at b = y2| J\- This level crossing 
is irrelevant to the low-T behavior, however, but modifies 
the M(B) curve somewhat at finite T. 



B. AFM si = 5/2 dimers with strong anisotropy 
interactions 



We now consider some cases of strong anisotropy in- 
teractions in the s\ — 5/2 dimer. In Fig. 22, we 
present low-T plots of M(B) with Jb = Jd = c|J|, 
where c = 0.1, 0.2, 0.3, 0.4, 0.45 and 0.5, as indicated. For 
the weak coupling case c = 0.1, the steps have uniform 
height. With increasing c, the curves shift monotonically 
to lower b values. Slight deviations in the uniformity 
of the step height are detectable in the first two steps 
for c = 0.2, but the nonuniformity in the step height is 
progressively more pronounced for c = 0.3 and 0.4, re- 
spectively. However, for c = 0.45, the individual steps 
have completely disappeared, and arc replaced by what 
appears to be a single, broad step centered at b/J « 0.5. 
At c = 0.5, M(B) rises to its maximum value with a very 
steep slope at b = 0, not showing any evidence for any 
level crossings. 

Next, we investigated the effects of large J a . In Fig. 23 
we plotted M/7 versus r )B/\ J\ for J a = c\J\, with c = 0.1 
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FIG. 23: Plots of M/ 7 versus jB/\J\ at fc s T/| J| = 0.01 with 
J a = c\J\ for the AFM si = 5/2 dimer. The cases c = 0.1 
(solid), c = 0.2 (dashed), c = 0.3 (dotted), c = 0.4 (dot- 
dashed), c = 0.5 (short dashed), and c = 0.6 (short dotted) 
are shown. 
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FIG. 25: Plots of C v /k B versus jB/\J\ at k B T/\J\ = 0.01 
for the AFM si = 5/2 dimer. The cases Jb = Jd = c\J\ 
for c = 0.45 (solid) and c = 0.5 (dot-dashed), J a = c\J\ for 
c = 0.4 (solid) and 0.5 (short-dashed) and J c = 0.5 1 J| (thick 
dashed) are shown. 




FIG. 24: Plots of M/j versus jB/\J\ at k B T/\J\ = 0.01 with 
J c = c\J\ for the AFM si = 5/2 dimer. The cases c = 0.1 
(solid), c = 0.2 (dashed), c = 0.3 (dotted), c = 0.4 (dot- 
dashed), and c = 0.5 (short dotted) are shown. 



(solid), 0.2 (dashed), 0.3 (dotted), 0.4 (dot-dashed), 0.5 
(short dashed), and 0.6 (short dotted). We see the the 
effects of strong J a are very different from those of strong 
Jb and Jd shown in Fig. 22. Instead of the step positions 
decreasing monotonically with increasing J b = J dl as J a 
increases, the highest three step positions decrease with 
increasing J a , but the lowest two step positions increase 
with increasing J a . 

We now consider the case of large J c . In Fig. 24, 
we plotted M/7 versus b/\J\ at k B T/\J\ = 0.01 for 
the si — 5/2 AFM dimer with J c = c\J\, where c = 



0.1, 0.2, 0.3, 0.4, and 0.5. These results are shown respec- 
tively as the solid, dashed, dotted, dot-dashed, and short 
dotted curves. Unlike the strong global anisotropy case 
pictured in Fig. 22 and the strong local axial anisotropy 
case pictured in Fig. 23, these curves do not shift signif- 
icantly to lower b values with increasing J C /\J\, but the 
step shapes are greatly altered. With c = 0.4, the first 
two steps are hard to discern, but rounded remnants of 
the three last steps are evident. At c = 0.5, the third step 
is hard to detect, and the remnants of the last two steps 
are reduced in magnitude. We remark that the curve 
with c = —0.3 is indistinguishable from the dotted curve 
for c = 0.3, so that in this case, strong J c /J of either 
sign can grossly alter the M(B) step behavior. 

In Fig. 25, we show Cy/fcs versus b/\J\ low-T curves 
for some of the AFM si = 5/2 cases pictured in Fig. 
22-24. The solid and dotted curves on the left-hand 
side of the figure correspond to Jb = Jd = c\J\ with 
c = 0.45 and 0.5, respectively, the solid and dot-dashed 
curves in the central portion of the figure correspond to 
Ja/\J\ = 0-4 and 0.5, respectively and the dashed curve 
running throughout the domain pictured corresponds to 
J c = 0.5 1 J\. We note that for the global anisotropy 
case shown, strong anisotropy pushes the Cv{B) peaks 
to lower b, and squeezes them together, so that they 
tend to overlap. At the limiting case Jb = Jd = 0.5 1 J|, 
the individual level crossings have been eliminated, and 
all ten of the peaks are combined into a single Schot- 
tky anomaly. For the local axial anisotropy interaction, 
strong FM anisotropic interaction squeeze the Cy(B) to- 
gether in the middle of the domain, with the limiting 
case of a single peak at J a /\J\ = 0.6. For the strong az- 
imuthal anisotropic exchange interaction case, however, 
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strong FM anisotropy shifts the double peak positions 
only slightly, and broadens the double peaks, so that only 
the highest set is clearly discernable, even at the very low 
T value plotted. In each case, the deviation of the peak 
height from the two-level prediction , Eq. I|41|). is sub- 
stantial, indicating more than two levels are important 
at low T, as for the three-level system with si = 1/2 
and Jb = \J\. Thus, the different types of strong spin 
anisotropy lead to clearly different low-T Cv{B) behav- 
iors for si > 1/2 dimers. 

IX. SUMMARY AND CONCLUSIONS 

In summary, we solved for the low-temperature mag- 
netization and specific heat of equal spin s\ antiferro- 
magnetic dimer single molecule magnets, including the 
most general set of anisotropic spin exchange interac- 
tions quadratic in the spin operators. The magnetiza- 
tion and specific heat exhibit steps and zeroes, respec- 
tively, at the non-universal level-crossing induction values 
B 1 g Si (9,(j)), but the magnetization steps and their mid- 
point slopes, plus the two peaks surrounding the specific 
heat zeroes all exhibit universal behavior at sufficiently 
low temperatures to first order in the anisotropy inter- 
action strengths. Strong anisotropy interactions gener- 
ally lead to highly non-universal behavior in antiferro- 
magnetic dimers. Local (or single-ion) anisotropy in- 
teractions lead to low-temperature magnetization step 
plateaus that have a much richer variation with the mag- 
netic induction B than do those obtained from global 
anisotropy interactions, provided that s\ > 1/2. For the 
most general quadratic anisotropic spin interactions at an 
arbitrary B, we derived simple, asymptotic analytic ex- 
pressions for the low-temperature magnetization, specific 
heat, inelastic neutron scattering cross-section, and elec- 
tron paramagnetic response susceptibility, which are ac- 
curate for weak anisotropy. We also derived an accurate 
expression for the level-crossing induction B\ c Si (9, </>), en- 
abling fast and accurate fits to experimental data. 

There were two low-T M (B) studies of Fe2 dimers. [lil 
[rH | For /i-oxalatotetrakis(acetylacetonato)Fe2, all five 
peaks in dM/dH were measured in pulsed mag- 
netic fields H. These evenly spaced peaks indi- 
cated little, if any, spin anisotropy effects. [lq On the 
other hand, studies of the first 2-3 dM/dH peaks 
in powdered samples of [Fe(salen)Cl]2, where salen is 
N, iV'-ethylenebis(salicylideneiminato), were much more 
interesting. 0| These data showed a broad first peak at 
B = 17 — 20T that was only partially resolvable into 
two separate peaks, followed by a sharp second peak 
at B = 36 T, consistent with local axial anisotropy of 
strength \J a /J\ as 0.1, as obtained from the derivatives 
of the curves shown in Fig. 13. From Eqs. ljfi!3 |l -l|68 jl . 
one could also have \(J a ± 3J e )/ J\ as 0.1. These amounts 
might perhaps be combined with a smaller \J C /J\, ob- 



tained from the derivatives of the curves pictured in Fig. 
14. 

Without a detailed single crystal study with magnetic 
fields along different crystal directions, it is impossible 
to determine the relative amounts of J a and J e that 
would best fit the data. However, the existing data on 
[Fe(salen)Cl]2 appear to be inconsistent with a predom- 
inant global anisotropy interaction of either type, as ob- 
tained from the derivatives of the curves shown in Figs. 
11 and 12. Only single crystal studies could determine 
if a small amount of such global anisotropy interactions 
were present in addition to one or more presumably larger 
local spin anisotropy interactions. In comparing the two 
materials cited above, it appears that the interaction of 
a Cl~ ion neighboring each Fe 3+ ion leads to strong local 
(or single-ion) anisotropy effects. In order to verify this 
hypothesis and to elucidate the details of the interactions, 
further experiments using single crystals in different field 
orientations on this and related Fe? dimers with 1-3 simi- 
larly bonded CI - ions are urg ed.ha^Oj We also urge sin- 
gle crystal data on some of the si = 9/2 [Mn,^ dimers, 
[LJ, |2lj, |22j , as well as on si = 1/2 dimers lacking in pre- 
dicted local spin anisotropy effects. 0, 0, 0, H4]To 
aid in the fits, we derived simple, useful formulas for 
the magnetization and specific heat at low temperature 
and sufficiently large magnetic induction. More impor- 
tant, we derived accurate analytic formulas for the elec- 
tron paramagnetic resonance susceptibility x(\bmB,ui) 
and the inelastic neutron scattering S(B, q, u), which al- 
low for a precise determination of the various microscopic 
anisotropy energies. 

With local anisotropy interactions, the total spin s is 
not a good quantum number, potentially modifying our 
understanding of quantum tunneling processes in single 
molecule magnets. It might also be possible to fit a vari- 
ety of experimental results using a smaller, consistent set 
of model parameters. |6j We emphasize that the study of 
dimer single molecule magnets, for which the most gen- 
eral anisotropic quadratic exchange interactions can be 
solved exactly, may be our best hope for attaining a more 
fundamental understanding of the underlying physics of 
more general single molecule magnets. 
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APPENDIX A 
Single-ion operators for arbitrary Si,s% 

For completeness, we present the results of the single- 
ion operators for general dimers consisting of magnetic 
ions with spins s\ and s 2 . We find 

+DT + x, SUS2 W+i)), (H7) 

--(-^(^..jcr) 

-a+ixjcjr>), (us) 

C™ llS2 = %,. 1)aa >/(a-m)(s-m-l), (119) 
^™ liS2 = //,.. \ - iu-. (120) 



'[(,s 1+S2 + l) 2 - S 2 ][ S 2 -( Sl - S2 ) 



y (4s 2 - i)s 2 

gl(gl + 1) - g 2 (s2 + 1) 

s(s + l) 



(121) 
(122) 



Note that for si = s 2 , £ s . Sl . Sl = and i] s , Sl , Sl reduces 
to Eq. (|22|) for ?j S)Sl . Also, note that these relations 
preserve the global spin properties of S z = S\ >z + 5 2 . z 
and S ff = Si )tr + S 2 ,a, S z \ip?) = m|Vr) and S a \^f) = 
A^ m \ip" l+,T ), as for equal spin dimers. In addition, the 
spins on asymmetric dimers with si 7^ s 2 would be ex- 
pected to also have different local anisotropy interactions, 
so that TL a and TL e would become 



K 

H' 



El S* 2 

i=l 
2 

^ ] Je,i(S i x — S i y ). 
i=l 



(123) 
(124) 



Specific heat details for s\ = 1/2 

We first present the numerators of the exact expres- 
sions for the specific heat with si = 1/2 and B\\i for 
i = x,y, z. We have 



1 



P J x ,y 



+F Xi ysmh(fiF XjV )e 
x [(J + J x , y )e~ 0J + (J x , y - 2 J B ,x)e a/SJ »'"] 



cosh(/?.F 3 



x([(J + ^,,) 2 + ^ 2 

A4 = F 2 cosh(2/3F 2 ) + ^e- /3 ( J+2Ji >) 

+ icosh(/3 J F z )(A,(J 2 +F 2 ) 
+e -/3(j+j b ) J(J + 2Jb ^ 

+F Z sinh(/3F 2 ) ( J b A, + J e -^ J+Jb > 



(125) 



(126) 



Eigenvalues for si = 1 in the crystal representation 

In the remainder of this appendix, we provide some 
details of our exact results for si = 1. The cubic equation 
for the three s = 1 eigenvalues is given by 

e n = -J-J b -J a + A„, for n = 2, 3, 4, (127) 
= -A n -(J Q - J b )\l + X n [b 2 + {J d - J e f] 
+ (Ja-Jb)[b 2 cos 2 + (J d - J e ) 2 } 
-{Jd - J e )(6 2 sin 2 (9cos(20). (128) 

For B||z, the cubic equation is easily solved to yield 



\n = Jb-Ja,±Vb 2 + (Jd-Je) 2 - (129) 

For -B||x, y, we have 



—2 



A„ = ±(J e - J d ),-J ytX ±^b 2 + J xy , (130) 



(131) 



where (J y ) corresponds to the upper (lower) sign. 
The six eigenvalues for the mixed s = 0, 2 states satisfy 



e rl — — 3J — — J a — 2Jb + A Tl 



We first define 



a = — J a 



d 3 = \l-d 



b 3 = s]-b±, 

d = -V6(J d + J e /3), 

e -Aj 

- _ 2 



(132) 

(133) 
(134) 

(135) 
(136) 
(137) 
(138) 
(139) 



and 
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QP = n(J b + J a /3)+pbcos9. 



(140) 
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Then the A n are the eigenvalues of the Hermitian matrix 
M given by 



M 



(Qzl 

d 





V 








d 

Q° 
K 
d 



o 

d 3 
h 

Qi 








d 

b x 

Q 2 -2 

e 



e 

a 


e 



<2° 



.(141) 



3 J J 



The resulting sixth order polynomial for the A„ is given 
elsewhere. [25) We note that for B\\i, M is block diago- 
nal, breaking up into matrices of ranks two and four. 
These cases are discussed in detail elsewhere, jl^ Here 
we only present the simplest data for which the eigen- 
values are determined either by linear or by quadratic 
equations. 

Simple special cases 

When only one of the Jj 7^ 0, the eigenvalues for B\\z 
simplify considerably. For J b 7^ 0, we have 

\ 7 n = 3 J + 2J b , 2J b , J b ± b, —2J b ± 2b. (142) 

For J c i 7^ 0, we find 



Ki = n,3J,±,JP + 9J*,±2^P + 3J*. (143) 
For J a 7^ 0, the eigenvalues can also be found analytically, 



K = - 2 ^±2b, J f±b, J f + ^±^j\j* + J*-JJ a . 

(144) 

We note that the ground state energy in this case is 



El = -\j-Ja-\J\j 2 + J 2 a -JJa, (145) 

which explicitly involves mixing of the s = m = and 
the s — 2, to = states. 

At the first level crossing with B\\z, we have for J a 7^ 
and the other Jj — 0, 

jB\ c f = I(j+[9J 2 +4J a 2 -4JJ a ] 1 /2). (146) 

We note that the first level crossing for J b ,Jd 7^ for 
B\\i is equivalent to that of si = 1/2, given by Eq. (|57|1 . 
At the second level crossing, simple formulas are only 
obtained for B\\z with one Jj 7^ 0. For B\\z and J a 7^ 0, 
J b ^ and Jd 7^ 0, respectively, we have 



t->1c,JS 

1 B 2 1 



-2J-J a , 

-2J-3J b , 



i (20 J 2 - 30 J 2 + 8[4J 4 - 6 J 2 J 2 ] 1 / 2 



(147) 

1/2 



Rotation to the induction representation 

The rotation from the crystal representation to the in- 
duction representation is obtained from 





(148) 

leading to B — Bz' . 28] This operation is equivalent to a 
rotation by — 7r/2 about the z axis, a rotation by 8 about 
the transformed x axis, and then a rotation by tt/2 — <fi 
about the transformed z axis.|2ij In effect, in using the 
above rotation matrix, we made the arbitrary choice that 
the rotated z axis lies in the x'z' plane. After the above 
rotation, it is still possible to rotate the crystal by an ar- 
bitrary angle x about the z' axis, keeping -B||i'. Hence, 
there are in effect an infinite number of equivalent rota- 
tions leading to B = Bz' . The resulting Hamiltonian 
matrix will then have off-diagonal elements that depend 
upon x- However, all such rotations necessarily lead to 
the identical, ^-independent, set of eigenvalues of the re- 
sulting diagonalized Hamiltonian matrix. We have ex- 
plicitly checked that the above rotation gives the exact 
cubic expression, Eq. (|4*7jl . for the si = 1/2 eigenvalues, 
and also leads to the correct eigenstate energies second 
order in each of the Jj for si = 1. We also showed explic- 
itly that x does not enter the eigenstate energies second 
order in J b for arbitrary s, Si,m. 



Global Hamiltonian 



The global axial and azimuthal anisotropy interactions 
in the rotated frame are 



-J b (s 2 z , cos 2 6 + S 2 , sin 2 (9 
-sm(2e){S x ,,S z ,}/2), 



(149) 



H'd - -Jd 



cos(2^)(5 2 , 



cos 



S 2 , sin 2 9 - S 2 , 



+ sm(20){S x /,S z /}/2 
- sin(20) (cos 6{S X > , S y > } + sin 6{ S y < , S z . } 

where {^4, B} = AB + BA is the anticommutator. 



(150) 
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Global Hamiltonian matrix elements 



and 



The operations of the rotated global anisotropy inter- 
actions upon these states may be written as 

KW?) = -f [(4m 2 + 2[ S ( S + l)-3m 2 ]sin 2 fl)|^) 

-sin(20) Y, (2m + a)A™\tp™ + °) 



o=±l 



-sin 2 9 Y F! m W. 



m-\-2cr\ 



(151) 



and 



(2sin 2 6icos(20)[3m 2 - s(s + 1)]|</?" 1 ) 



-2 sin 9 E (2m + a)A° s 



<r=±l 



x [cos 9 cos(20) + iff sin(20)] \ip™ +a ) 

— pom ^ + CQS 2 ^ cos ( 2 0) 



cr=±l 



+2icrcos6»sin(2 ( /))]|^™ +2 ' T ; 



(152) 



where is defined by Eq. (|17|) . We note that in this 
representation, both Hb and Hd preserve the global spin 
quantum number s, but allow Am = ±1, ±2 transitions. 



Local Hamiltonian 

With regard to the local spin anisotropy terms in the 
rotated coordinate system, we write 

H' a = -J a (Oi cos 2 9 + 2 sin 2 9 - 3 )l53) 



and 



H„ = —J. 



cos(20) sin 2 9 + e> 2 cos 2 

+ isin(20)O3-O 4 ) 
sin(20) (o 5 cos 9 + O e sin ( 



where 



i=l 
2 

i=i 

2 

i=l 
2 

i=l 
2 



(154) 

(155) 
(156) 
(157) 
(158) 
(159) 



@6 — ^ ^ {Sjy'Sjz' ^rSiz'Sjy'). 



(160) 



Local Hamiltonian matrix element components 



The operations of these interactions are given by 



+ e #;xV/w>), (i6i) 

a'=±l 



Er crm 



a' = ±l 
m+2a\ 



a=±l 



+ E K Z/\^))^ (162) 

<T=±1 

- e ^'^ ,CT Vr + + 2^)), (163) 

cr'=±l 



'■p 



a' = ±l 
m+2a\ 



a=±l 



<T,rr' = ±l 

o=±l 

+ E ^ ct V s + + 2 ^>), (165) 

cr'=±l 

= Ti E (^>r +CT > 

CT = ±1 

- E ^rWw)), (we) 



er'=±l 



where 

M s m si = -Va s , Sl +4[ S ( S + l)-l](l-a s , Sl ), 



(167) 



\jm,o _ m-(cr +l)/2 r1 o-cr m+(o- -l)/2/-. rq n 

JV s,s'i - 2^ U s+(<t'+1)/2,si U s+(3 CT '+l)/2,si ^! iD8 J 
cr=±l 



P-f = 2Ar(2m + ff)a SiS1 , 
and 

nm,a,a' _ ft— moo' — (</+l)/2 j^m+o 

rt s,si — °s+(cr' + l)/2,si s+(3<r' + l)/2,si 



(169) 



r ,-mo-cr'-(cr' + l)/2 nm M 7rA 

+ °s+(3 C r' + l)/2, Sl iy 8 +(ff'+l)/2,s I > U' u ; 
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where r? s , si , G™ S1 , i^ 7 , K*£ , and L* S1 are given by 
Eqs. (123) and l|2l ^ -l|28 jl . respectively. We note that for 



o' = ±1, N?£ = 2H%° 



APPENDIX C 
Inelastic neutron scattering cross-section 

Here we outline the derivation and provide the im- 
portant results for the low-T inelastic neutron scatter- 
ing cross-section, S(B,q,u). At low T, we set E™ Si w 



E. 



m,(0) 



E™s{ 1] given by Eqs. and (J7TJ- Since this 
energy was evaluated in the induction representation, it is 
easiest to evaluate S(q,u>) by choosing the axes e a = i', 
where i — x,y, z. The i' can be obtained from the i using 
the inverse of the rotation matrix in Eq. |148jl . Then, 
for q = (sin# 9 cos <fi q , sui9 q sin <p q , cos9 q ) in the crystal 
representation, we have 

q = (sin 6b, q cos cf> b>q , sin 6 b<q sin <f> b<q , cos 9 biQ ) (171) 

in the induction representation, where 



cos ti b>q 
sin 9 b<q cos 4> b , q 
sin06 >g sin 4> btq 



sm y sm t) cos 



sm t/„ cos V cos 



g \y q f , 

g ) — sin # cos , 
(172) 



Then, expanding the wave functions to first order in the 
anisotropy interactions, 



Irs 



\9» 



|^')((^'|iVW 



(173) 



' cr=±l 



Eyym,<r 

cr' = ±l 



_ym,<7| m+2en 
2 s.si Irs / 



J[o-'(2s + 1) + 2] 



I , ^TTi \ 



■\sm.(T,<T 



^\ajB + J[a'(2s+l) + 2] 
2ct 7 B + J[a'(2s + 1) + 2] '^+ 2ct ' V 



(174) 



where 



i(2m + a)Ar[sin(20)(j b S ! f 
-J d % sl cos(2<^ 



-2ioJ s ^ sin#sin(20) 



(175) 



_ TTiam 



y 6.a ' 



+ ^,e 1 ( 1 + COs2 COS(20) 

+2iaJ^ 1 cos#sin(20) , 



(176) 



i 



and 



~-H™£ [J a cos 2 6» + J e sin 2 9 cos(2</>)] 

+ijV£'f ' sin 2 9[J a - J e cos(2$], (177) 

^Kf/ (sin(20)[J Q - J e cos(20)] 
-2ia\7 e sin0sin(2<£)Y (178) 



J a sin 9 



+J e ((l + cos 2 0) cos(20) 



+2icrcos0sin(2</>) 



(179) 



where fl^af' an d ^f.'fi' are given by Eqs. |2"6"fl and QZJj 
respectively, iV^'f' and Rff^' are given by Eqs. fTH^I 
and l|170p . respectively, and J^ 1 and Jj'^ 1 are given by 
Eqs. O and ffify. respectively We note that for ct' = 

±l,Ar™.-' = 2^'. 

Now, we evaluate the matrix elements 



*C«(<l) = (^|5t(q,0)|^), 



(180) 



including the leading corrections to the wave functions 
due to the anisotropy interactions. We then find to first 
order in the anisotropy interactions, 



isin(q-d) 



<j'=±i 



X (D s+ ( ct / + l)/2, Sl ^ 's ,s+a' An,m 

+D'T Hal+l)/2 , s M%\^,) 

+Df +{a , +1 y 2 ^ +a ,\N£l)), 



(181) 

K£,(q) = \Y, Vos{q-d)(Al m &- s , s 5rn >m+a 

<r=±l - 

— icr sin(q ■ d) (— er') 

cr'=±l 

( r -aa'rn-(a' + l)/2 r_ r_ 

y s+((T / + l)/2,si (, s,s+d'"m,in+i7 



r ,-crcr'm-(<j' + l)/2 , , .(1) | m+a\ 
+ °s+(<j' + l)/2, Sl \^-s.Jn\^s+a'l 

_ r ,aa'm-(a'+l)/2, m-tr\ *r(l) \ 
°s+(<t' + 1)/2, Si \rj+a> | JV s,m/ 



(182) 
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and (q) is obtained from M™J.™,(q) by multiplying 

the entire quantity inside the summation over a by —ia. 

Now, it is convenient to break up the contributions 
from |Mj^?/(<z)| and \M™'™, (q)| into their symmetric 
and antisymmetric parts, with 

|2 



M™' m (q) 

I s,s, sym \^ / 



1 / , M m,m , x|2 | M _m,m ( y2 



(183) 



etc. Then, the combined contributions from l-M^™'™ (q)| 2 

and |-^™ s '™ ym (<z)| 2 arising from the first order perturba- 
tions of the wave functions in addition to the Hartree 
approximation are 



S[ U \B t q,u>) 



2si s 

£ £ e_ 

s—0 in— — s 



x cos 2 (q-d) sin 2 O btq JT£(w, 9, 



■ sm 2 (q ■ d) sin 2 9 b;g J 7 ^f i ((J, 9, t 



■ cos 2 (q • d)'- 



sm 



rr-ra.s { 



where we suppress the oj, 9, 4> arguments of the J r ™fi^ 
functions for simplicity of presentation. The contribu- 
tions Tn^f^^uj , 9 , (f)) second order in the anisotropy in- 
teractions may be written as 



/7-m 
•*~ 1,. 



S (2) _ 



E 

<T,<7 1 — ±1 p— p'=0 

ct,(t'=±1 p=0 p'=0 



^; (2) = e ee^kc-cu 



_ pm+pa n 

8, SI S + 2pV,Sly 



(191) 



S'l s+(2p' + l)er',si/ 



X„4 : 



2 V^Wi 

3 1 



(192) 



a.tj'— ±1 p— p'— 



x„4 : 



3 1 



(193) 



^";; (2) = E EE 5 K^-« +1K J 



cr.cr'— ±1 p— p'— 

X.A?' P '(M 



(194) 



+ sin 2 (qr • d) 



2 - sin 2 9, 



b,q 



(184) 



where 



£j™ W (w,M) (185) 



are the contributions to the ^-" n through second order in 
the anisotropy interactions, respectively. The terms to 
zeroth order in the anisotropy interactions are given by 
Eqs. I|82l) - (|85|l in the text. The non- vanishing first order 
terms are 



T. 



2,si 



ft 



J2 6(u; + E™ 1 -E i 



s+a' ,sij 



cr'=±l 



xD 



s+(a' + l)/2, Sl Q2 i 



T 



■m,s(l) 

4,81 



Q2 



Q7 



= -n 



(186) 



E S (" + E Z 



pm+a > 
-^s+a'^i/ 



E 



cr,CT'=±l 

A-t— (f'+l)/2— crcr'm^o-,0-' 
XO s+(ct' + 1)/2,s 1 ^4 



(187) 



E 



rrm,<7<7 

s+(j' + 1+2ou')/2,si s+J , (l-J)/2,si 

2[ct'(2s + 1) + 1 + ct] 

f^<7<7 m+(er — 1)/2 



(188) 



a=±l 



f70 s+(o-' + l+2ao'')/2,si s+ct'(1-ct)/2,si 

2[er'(2s + l) + l + cr] 



(189) 



(!)/ 



J a - 3[J a cos 6» + J e sin 2 9 cos(20)] , 



(190) 



We note that we have suppressed the a, a' ,m, s, s% de- 
pendencies of the A v ^ v (9, (f>) for brevity. These functions 
are given by 



A° { ° = A*' 1 = 0, 



A{'° = [tf\9,<P) + 2a StS JV 



2 V^: ' 



+< S1 / 3 (2) (M)] (GUY, 



A 2 {° = [f^ ) (9,cj ) ) + 2a^ 1 fP(9, l 



A 



+<8 1 /f (M)] (G! 
= /f (M)Kf) 2 , 

^' 2 = /f (MK^) 3 , 

n.i!!i)/2 _ ^(2)(e )T 4)/^±ff',<r'^2 
f(2)/ 



,2) 2 , 



(195) 

(196) 

(197) 
(198) 
(199) 



(200) 



4'° = (<£;?) *lf? ) {e,4>) + 2a s , s j, 

,2 f( 2 ) 



+<8 1 /r j (^^)] + (G 2 ff ;r- <T )7 3 v- 



;Z'- ff ') 2 /3 (2) ( 

(2) //a , „, /-(2), 



+2G2J G%' t Z ' a [f 2 *'{0, <fi) + a StSl f^ yw, V)i , 

(201) 

A'" = (G^') 2 [/i 2) (0,0)+2a_ 1 / 5 (2) (^, r i) 

+<sA 2) M+ipli'^'ff^ (M) 
+2G^'G^' [/f } (0, 0) + a S)Sl /f 0)] , 

(202) 

= / 3 (2) (0,0)(G^'^') 2 , (203) 

(204) 



2 

2.1 



A 



A 



1±1,0 



(2), 



+ 2^ Sl / 2 (2) (0,, 
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[fi 2 \e^) + 2a SiSl fP(9, ( p) 

+< si /f (e^Wlt) 2 ' ( 206 ) 

/ 3 (2) (^0)(G^'' ±CT ) 2 , (207) 

(Gl%) 2 tf\e^) + (Gl%>-°) 2 f (2 \e l( j>) 

+2G^'G^'-V 9 (2) (e,0)], (208) 



(205) 



A 



2±1,0 



«1±1,1 



. 1,1 



.3,1 



(209) 



+2a s , Sl f¥\9^)+ a i Si f^(e, ( j ) )] 

I c\^1<7,(7 ,i<7 fi(T,(J ,ifT, — <7 

' Z(j "4.1 ^4.4 



x[/f ) (0,0)+« SjSl /f , (^ 



(2), 



(210) 



a^ = (gi%-°Y[l 



^ r 1 2) (^ 



+{Glf-^- a 'ffi 2 \e^) 

I <r\f^l<7,<7 , — <7 ^r<7,<7 , — (7, — <7 



x[/ 6 (2) (M) +«.,.! / 6 W (* 



(2)/ 



+ (G^'- CT )7r J (^0) + 2G^- CT 
x (g!; 2 ct ' ct [/ 8 (2) (0,0) + a a , ai / 9 (2) (0, 0)] 



+G££''-*'-*7< 2) (M)), 



(211) 



A 3 4 ° = (G^' ff ) 2 [/i 2) (0,0) 

+2a S;S jf ) (^)+a 2 iS jf ) (0,0)] 



cr,o-',cr,-(7'^2 f (2) ^ 



( G 4,5 )7, 



r)X^Cr,Cr .(7^(7,(7 ,(7, — (7 



x [/ 6 W (M) +«.,.! /r(M)], 

4 ±M = /i 2) (M)(G£l' ,w ) 2 , 



(212) 
(213) 



4' 1 = (g:; 5 ct ''- ct ' ct ') 2 /6 (2) (^^) 

+ (G 4 ^'<< T ') 2 /f ) (0,0) 

+2G^' CT^' f^ 2) {6, 0), (214) 
= /f (M)(GI£ V,,7 ') 2 > (215) 
where the second order angular functions are given by 
sin 2 9 



f[ 2 \o 



/2 ([J 6 -cos(20)J d ] 2 cos 2 

+J 2 sin 2 (20)), (216) 

-j2- ([ J t ~ J d cos(20)][J Q - J e cos(20)] 



jf (M) 



sin 2 (9 



x cos 2 6» + J d J e sin 2 (20)]), (217) 

, , ([J a -cos(20)J e ] 2 cos 2 

+J 2 sin 2 (20)), (218) 
/fW) = ^(j b sin 2 + J d (l + cos 2 0)cos(20)) 2 

(219) 



-4^ cos 2 9 sin 2 (20), 
1 



fPi^A) = -^ 2 -[Jhsin 2 + J d (l + cos 2 0)cos(20)] 
x [J a sin 2 9 + J e (l + cos 2 0) cos(20)] 



cos 2 sin 2 (20), 



(220) 



/6 2) (M) = j2-(jasin 2 + J e (l + cos 2 0)cos(20))' 



-4^| cos 2 sin 2 (20), 



(221) 



/t 2) (^0) = j2-(ja-3[J a cos 2 + J e sin 2 0cos(20)]) 2 , 
f (2) {9,cp) = -l(j a -3[J a cos 2 + J e sin 2 0cos(20)] 



x J b sin 2 9 + J d {\ + cos 2 9) cos(20)j ,(222) 
/g 2) (^0) = j 2 -(ja-3[J a cos 2 + J e sin 2 0cos(20)]) 

x (j a sin 2 6> + J e (l + cos 2 (9) cos(20)) ,(223) 
and the constants are given by 



= J(2m + a) 1CTTO 



Gl,2 



Gcr,<x 
1,4 



G(7 ,(7 
1,5 - 



GCr,(T 

2,1 



2jB 

jpam 
4 7 S ' 



TD7n,(J.<7 

s, Si 



4(ct 7 B + J[<t'(2s + 1) + 2]) ' 



S,Sl 



4^2cr 7 B + J[cr'(2s + l) + 2]) 



Ja(2m + cr) 



E H™^- 



" ^ CT=±1 



n m+ CT (CT+l)/2 \ 
Xi ^ S +(l+ ( T')/2,siy' 



■oCT.ct' _ \ ~> n m+o-(CT+l) 

2 ' 2 ~~ 8-7B ^ CTi VKl+<7')/2,si 

' CT=±1 



A - r s-cr'(5 : -l)/2,si' 



(224) 

(225) 



(226) 



(227) 



(228) 



(229) 



2,3 



E 

CT=±1 



aD 



_ s+<j'(o r +l) + (a+l)/2,s 1 

4[cr'(2s + 1) + 2 + aa'(l- a)} 
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x H m a 



(230) G] 



1 ( ,— n rn+a(a+l)/2 

A 2^ \ <J - Ly s+(7'(a+l)+(<T"+l)/2 lSl 



4 \ 



i? 



s+o-"(1-ct)/2,si 



ajB / J + a'(2s + I + a" - a"a) + 2 J' 

(231) 



I " s +o-"(1-ct)/2,s 
ct=±1 V 



T-.m+<r(cr+l) 
s+(ct+1)<7' + (<t"+1)/2,si 



2ajB/J + a' (2s + 1 + cr" - o"o) + 2 J' 

(232) 

|^((2m + a')^ m +^' m 

-(2m + 2cr + crVs' (m+ ' T) ^s m ), (233) 

__^_( Aa'(m+2cr) perm 
8"fB V s 

_ i4 j'm f «r(m+<T')j ) ( 2 34) 

4 [a'(2 S + 1) + 2] ' (2d&J 

(7 — ±1 

jftiA^tryfl + J[a'{2s + 1) + 2]) 
xK?-< r "&- 1 V 2 w' ), (236) 



s + (cr + l)CT' 



- e 

F=±1 V8(2(77B + J[ct'(2s + 1) + 2]) 
x^ s T 1 +(1 ^ ) ' T<T " /2, ' T '^ (237) 

trV'J 



^[2™ + cr" + cr(l - cr)] 



2 7 B < 

r ,-crcr'm-crcr'cr"(o r +l)/2-(cr' + l)/2 
X ° S + (<T' + l)/2, Sl 

4 o-"m+ C r C r"(l-ff)/2 
X/1 s+o-'(l-ff)/2 ' 

(7 (7 J 



4,4 



- -i E 



/ /// pm+<r"(l — a)/2,a,a' 



aIX s+a"'(l-a)/2.s 1 

-=±1 ' 

cr"(7"'m-crcr"cr"'(CT+l)/2-(o-"' + l)/2 



c 



s+(<7"' + 1)/2+<t'(ct+1),si 



' ajB/J + a'(2s + 1 + cr'" - + 2/ ' 

(241) 



n o,a',o",a"' 1 I '"— T ^m+aa"(l-a)/2,cr' 

^4,5 - -o 2^ l CT <JjA S + CT "'(l-a)/2, Sl 



CT=±1 



C* 



Trm.a' r>m+a, — <7,—<7 
-"«,si ' rl s+2cr',s 1 



// /// / /// , * \ /n " III ! — i i \ 

-a a m—(a + l)/2 — <xct ct (<t+1) 
+ (<t"' + 1) /2jV (ff+1) ,si 

" 2a^/B / J + a' (2s + 1 + cr'" - cr"'cr) + 2, 

(242) 

where we have used the identities Az™ 1 ^ 1 = A™, 

rp—m — 2 T?m P/ m, ~ ° 

/c:;r'.a,:,i a- .V;;; . 

For completeness, we write the contributions to 
S(B,q,u>) from the off-diagonal and antisymmetric di- 
agonal matrix element terms that are first order in the 
anisotropy interactions. These could be measured by col- 
lecting data at specific (9b, q ,<i>b,q) values. We have 

2si s 

S 2 1C) = E E e" /3J! - 1 sin 2 ^ 9 cos(2^, g ) 

s—Q m— — s 

x (cos 2 (q • d) Y, + E Z, W) 

I <7=±1 

+ sin 2 (g-d) ^ 5 (u + E™ sl -E™+°J 



(243) 



(238) 



E_pCTm+crCT"(l-o 7 )/2 
CTt s+a'(l-a)/2 

^y— a'a"m— crcj' a" (a+1) — (a' + l)/2 
X ° S +(<t' + 1)/2 :Si 



' CT=±1 



2si s 

£ ^ e -^i S in(20 M )^cos 2 (<z-d) 

MM, 06,,)) 

x ^(w)G 6 ,i 

+ £ % + ^ ai -^+ <r ) G 6,2' 

CT=±1 

+ sin 2 (q • d) 



- E 



CT=±1 



„ H -m+ < r(l-cr)/2, < T' 
CT gj "s+g"(l-CT)/2, ; ii 

4[cr'(2s + l) + 2 + cr'cr"(l-CT) 



^-,-crcr"m-((j" + l)/2 
XL/ S + (cr" + l)/2 + (T'(CT + l),S 1 



(239) 



(240) 



]T 5(0; + ^-^^) 

'=±1 

x I (h^ (9, <j>, <j> b>q ) + a s>3l h { 2 1] (9, <p, cj> b ^ q ) 
xG( 3 + h 2 1 \9,4>,<p b , q )G( 4 } 



CT,cr'=±l 
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and 



x { (h^ (6, ^ ct> b , q ) + a a , Bl 4 1} {6, 0, 

(244) 

E E e _/3J! ^i sin 2 M sin(2&,,) 

s— m— — s 

x jcos 2 (g • d) X! <K W + ^" ! - ^0 

I <7=±1 

+ sin 2 (g-d) J] ^ + ^. 1 -^,,J 

CT,0-' = ±1 



(245) 



The angular functions and coefficients are given by 



/ 2 (1) (M) 

/ B (1) (M) 



J b sin 2 6> 



+ J d (l + cos 2 6) cos(20)j , (246) 
j(j a sin 2 # 

+J e (l + cos 2 (?) cos(2</>)) , (247) 

y cos<9sin(2</>), (248) 

Je 



cos 6* sin(2</>), 



(249) 



— (cos <fi b , g sin(20) [ J b - J d cos(2</>)] 
+2sin0 M J d sin6>sin(2</>)), (250) 
{0,4>,4>b,q) = j (cos <f> btq sin(26>) [ J a - J e cos(20)] 
+2sin0 6i gJ e sin0sin(20)), (251) 



J(2m + a)(A™f 
87S 



(252) 



Go, a 
5,2 



cJ „~oo'm-(l+o')/2 
"167,8 °s+(ff'+l)/2, Sl 

v \ ~> r ,oo'm+(o'+2oo'-l)/2 
X Cr ° S +(<r' + l)/2, Sl 



= ±1 



GO", (J 
5.3 



= ic 



pCrcr>n-(l-er)/2 
X s+<7'(1-ct)/2 ' 

T'm-(l+ff')/ 2 



(253) 



16 s+(^' + l)/2,si 



E 

CT=±1 



/ (T ^s+(l+o'+2oo')/2,s 1 



V2cr7B/J + cr'(2s + 1) + 1 + CT 



crcrm+(cr— l)/2,(Tcr / 
s-o'(o-l)/2,s 1 



^6,2 



^6.3 



^[s(s + l)-3m 2 ], 
a' J 



(254) 
(255) 
(256) 



D': 



8 7j g^s+(a-' + l)/2, Sl 

x [(2toct + <j),4 

CT,CT=±1 

-aa'm— (!+' 



<rcrm 

s+ct'(ct+1)/2 



XC 



s+(ct' + 1)/2 



^6,4 



1 g- Ly s+(CT' + l)/2,s 1 



7j<7 



-a<r')/2] 
si ] " 

E 



cr,CT=±l 
-(l-CTa-')/2 



(257) 



s+(1+(t'-2cto-')/2,si 



0-7.B + J[-ct'(2s + 1) + a - 1] 



xi? 



Tfl ,00 , — OO 



s+a'(a+l)/2, Sl 



(258) 



GO.O 
6,5 



'J 



87B 



c 



-crcr'm-(<7' + l)/2 
s+(cr' + l)/2,Si 



(2m + cr) 



E_ r) m+o-((7+l)/ 2 ,™ 
CT - L V|>' + l)/2, S iA.-a'(<F- 



l)/2. 



o=±l 



(259) 



GCT,cr 
6,6 



J 



C 



-oo'm-(o' + l)/2 



16 «+(<t' + 1)/2,si 



E 

CT = ±1 



— nm|ff(l-(r) / '2, oo. a o 
<jH s-(7'(-a-l)/2,s 1 ' 

ajB + J[a'(2s + l) + l + a] 



n m+<r(F+l)/2 
-^s+(1+<t'+2cto-')/2,si 



7,1 



GO", (J 
7,2 



47B 



(2m + ( r)( J 4r) 2 , 



(260) 



(261) 



C J s-i-(7<7'm-((7' + l)/2 
"^^ Lj s+(o'+1)/2,s 1 



o=±l 



—f~,Oo'm—(l — o' — 2oo')/2 

si 



\ ~> / — r <r<r m-(l-cr 



GO.O 
7,3 



pCTCTm-(l-o-)/2 
Xr s+o'(l-o)/2 J> 

1 ^-,-oVm-(cr' + l)/2 
g° s +( C r' + l)/2,s 1 



(262) 



I ct)/2, 



{l-o)/2.oo' 

si 



c 



-CTCT'm-(l+0-'+2CTCr')/ 2 N 

s+(1+o-'+2ctct')/2,si 

' 2ajB/J + a' '(2s + 1) + 1 + <r , 



(263) 
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Higher order corrections to the EPR response 

Here we evaluate the EPR response X-a a (-B,w) in- 
cluding the corrections to the induction representation 
wave functions first order in the anisotropy interactions. 
It can be shown that second order perturbations to the 
wave functions only contribute to the EPR response to 
third (and higher) order in the anisotropy interactions, 
and are hence neglected. From Eq. I|88|) . we expand the 
exact wave functions \<p n ) to first order in the anisotropy 
interactions. To this order, s and m are still good quan- 
tum numbers, so we rewrite \<fi n } — |</J™). We then find 



2*i 



2 s i 



s=0 m= — s s'=0 m' = — s' 

EX , X _ I 7~° a 1 2 

u s ,s u m %m+aa | | 

CT = ±1 

i s , \ x , I ya+a 

~r u s,s / u m' ,m+a+a' 

a'=±a 



' |2 



E ^s',s+2tr"|-^2<7"| 
a" — ±cr 



^ ] <5s',s+2cr" ( ^m',m+o-+£r' | Z^,7 
(7 1 ,(t" — ±er 



i r I 7 cr+2o-' |2\ 

i" m',m+tT+2tT' 1^20-" I I 



(264) 



where the amplitudes Z™_ s m are given to first order in 
the anisotropy energies by 



^ — a ® ( j^—2-\-amyrn,—a j\ am (^pn — a,a\ 



(265) 



(266) 



2 1 B\ s 



a am ( -\ ;m — a, a \ * \ 

A S [Vs.S! ) ), 



Aam i/A;m,(T A am ( ~\A j"i+(T, a 

A s+2a" VV sJ 1 - A s \ VV s+2a", Sl 



J[a"(2s + l) + 2] 

Z a + a ' — 

2<T " a' 1 B + J[a"{2s + l) + 2] 

( A a(m+a') v m,a' ,a" 
" 1 A s+2a" ^s.si 

Aaml ym+a+a' ,— a' a' 
\ A - s +2a".s 1 



^a + 2a' 

2a" 



(267) 



(268) 



(269) 



2a'jB + J[cr"(2,s + l) + 2] 



( A<?( m +2<?')^rm,a' ,a'' 



Aam /-i;m+cr+2o-',— a',— a" \ 
~ A s \-y s +2a", Sl ) 



(270) 



leading to 



Z, 



a+a' I 2 



J 



(27B )2(^ (m+ff ' , ( 2 ™ + -')^'" 

-A a s rn {2m + 2<j + <j')A a s '( m+ ^y 
x(/ 1 (2) (0,0)+2a^ 1 / 2 (2) (^,^) 

+< Sl /3 (2) (^' 



(271) 



Z n 



J\2m-af(A^) 2 (f{2) 



(4 7 5) 2 



(2), 



,™ ,2 = tf\e,4>) 

1 2CT "! 16[ct"(2s + 1) + 2] 2 



(/f (M) 



(272) 



CT=±1 



TT'm+a(\ — a) /2,a" \ 
s.si J 



I 2(7" I 



J 2 /f( 



|^'(T+2(T / I 2 

I 2a" I 



(273) 

) 

16(o-'7_B + J[cr"(2s + 1) + 2]) 

xf y aATtrl ( -^ )/2 i?r s +CT(1 ^ )/2 ' CTV ' ,N ) 2 , 

\ / ; S + CT"(l + <7) S,Sl I ' 

CT=±1 

(274) 

J 2 /f(M) 



64 



(2a' 7 B + J[ct"(2s + 1) + 2]) 

/ — Aam+aa' (l+a) j^a'm+aa' (l-a)/2,a" 

\ s+<t"(1+ct) ^s,si 



CT=±1 



(275) 



where £ s m Sl = £ s m,(0) + 5 is given by Eqs. J7U) and 
(ED, A™, F*, Kfg', a s , Sl , and R™f^' are given 

by Eqs. Q, JTZJ,'©, and l|T7Q|l . respectively, 

the fn\0, 4>) are given in Eqs. l(2TH|l - (|2^|l . and we have 
used the relations following Eq. Q242J1 ■ We note that since 
fj = ±1 for clockwise and counterclockwise oscillatory 
induction senses, setting a', a" = ±1 is equivalent to 
setting them equal to ±cr. 



For the three cases with 



to = m. 



there arc no asso- 



ciated resonant magnetic inductions, and since the EPR 
frequency u> cannot be easily varied in an experiment, 
these terms are irrelevant. For the 10 remaining cases 
with to' 7^ m, the resonant magnetic induction is then 
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found to have the form given by Eq. (|92() . The term 
with amplitude Z 2a has a\ = 1/2, b\ = 2(m + a), and 
ci = 0. The term with amplitude Z^ a has 02 = 1, 
b 2 = (2m — cr)/2, and c 2 = 0. The two terms with 
amplitudes Z% , f where a" = ±cr have 03 = 1, 63 = 
cr[2 + er"(2s + 1)], and c 3 = -(2m + a)/2. The two terms 
with amplitudes Z 2 °„, where a" = ±a, have 04 = 2, 
64 = 2cr [2 + cr"(2s + 1)], and c 4 = 4(m + a). The four 
terms with amplitudes Z^if" , where cr',cr" = ±er, have 
a 5 = 1/(1 + 2aa'), b 5 = [2 + a" (2s + l)]/(a + 2ct')], and 
C5 = to + a' + a/2. 



APPENDIX D 
Second order induction representation Hamiltonian 

In this appendix, we evaluate the corrections to the 
eigenstate energies second order in the four anisotropy 
interaction energies Jj for j = a, b, d, e. The operations 
of the rotated Hamilatonian Ti' upon the eigenstates \<p™) 
may be written as 



ay' I , ~ r < 



(7=±1 



/ 7 jm.cr I , „m+er\ . -\ jm.cr 1 , „m+2cr\ \ 

<**,<.! l ( / 3 s + 2 <W 

+^r' ff Vr + + 2 2 CT ')), (276) 



cr,er' = ±l 



where E™s{° ] and £™^ X) are given by Eqs. (J23J and (JHJ, 
respectively, and the 1A to y functions are given by Eqs. 
(|175[) to (|179[) . respectively. 



Second order eigenstate energies 

From Eq. (|276|) . the second order eigenstate energies 
may be written as 



E m,(2) = J_ J 



a=±l 



2 1 

+ 2 



S,S\ 



\w. 



m,cr' 12 

S,Si I 



j[2 + (2* + iy] 



^rn,o",fT |2 



7 ^ i VJ[2+(2s + 1) ( t']+ct 7j B 

m,cr,a | 2 . 

)• (277) 



I'-' S,Si 



J[2 + (2s + 1)ct'] + 2(77^ 
For simplicity, we rewrite this as 



nm,(2) _ ™,(2)M | nim,(2)V , pm,(2)W 



B m,(2)V 



£1 



where 



(278) 



m,(2)W 
s,si 



m sm 



[4s(s + 1) - 8m 2 - 1] 



2 1 B 

x(cos 2 0[J fc s ; Q Sl -cos(20)J rf s ;: 1 ] 



W(20)(J d \ Sl )- 



[2s(s + 1) - 2m 2 - 1] 



(279) 



87S 

x ([sin 2 9 J*'* 1 + (1 + cos 2 0) cos(20) J s d 



s,sii2 



pm,(2)W 



™,(2)^ = 



+4 cos 2 sin 2 (20) (J^ 1 ) 2 J, (280) 

d m / \ 2 

s,Sl ' J a - 3[J a cos 2 6» + J e sin 2 6>cos(20)] 



16J 



(281) 



fZJqB/J) sin 2 9 
2J 

-x([Ja- J e COs(2(f>)} 2 COS 2 9 



+ sin 2 (20) J\ 



(282) 



64J 



( Ja sin 2 
+J e (l + cos 2 9) cos(20) 
+4 cos 2 9 sin 2 (20) J 2 



(283) 



(a 2 - m 2 )[(s - l) 2 - m 2 }^ 2 ^ 2 ^^ 
(2s - 1) 

+ r ls+2 : s 1 7 ls+l,s 1 

[(s + l) 2 - TO 2 ] [(s + 2) 2 - TO 2 ] 
X (2s + 3) ' 

Vs^xVs — l^! ( g2 TO 2 ) 

(2s - l) 2 - X 2 
x((2s-1)[(s-1)(s-2)+to 2 ] 

— to(2s — 3)a;^ 

^+2, 5l ^+l, Sl [(^+l) 2 -^ 2 ] 

(2s + 3)2 - x 2 
x((2s + 3)[(s + 2)(s + 3) +to 2 ] 

-m(2s + 5)A (284) 

^jsiVs-^i ( (2s „ lU m 4 

(2s-l)2-4x 2 V^ S L)[m 
+to 2 (6s 2 - 18s + 11) 
+s(s-l)(«-2)(a-3)] 
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-Amx(2s - 3)(s 2 - 3s 

s+2,s 1 r 1s+l,s 1 



(285) 



(2s + 3) 2 - Ax 2 



(2s + 3)[m 4 



+to 2 (6s 2 + 30s + 35) 

+ (s + l)(s + 2)(s + 3)(s + 4)] 

-4ma;(2s + 5)(s 2 + 5s + 5 + m 2 



(286) 



There is a remarkable amount of symmetry in the an- 
gular dependence of the eigenstate energies. We note that 
i£™sP (0) 4>) an( i E™s[ 2S>U (0, 4>) have the same forms, dif- 
fering in the replacements of the interactions J^ 1 and 
Jjg 1 with J a and J e , respectively, and with different 
overall constant functions. The same comparison can 



also be made with E™s[ 2 ^( 



and £ S "M 2)V In 



addition, we note that there is a remarkable similarity 
in the 0, <fi dependence of E™s^ W with that of the local 
spin anisotropy part of B l £{\\0, </>) given by Eq. (|MJl . 

The contributions to the sth level crossing second order 
in the anisotropy interactions are calculated as indicated 
in Eq. lJ9fj)). and are found to be 



B = -Js/j 



jY,a n (s, Sl )fW(9,<f>), 



71=1 



(287) 



where the fn^ are given by Eqs. H216l) - (|222l) . respec- 
tively. 



Second order level crossing coefficients 

The coefficients are given by 
(9 - 20s + 12s 2 ) 



<n(s,si) 



2s 



(288) 



/ s Mj +4si(si + l)a 2 ,i(s)] , . 

tt 2 S,si = , : 1U0 : , (289 

s(2s + l)(2s + 3) 

02,0 («) = 3(9 + s-23s 2 + 4s 3 + 12s 4 ), 

a 2) i(s) = 9 -8s -4s 2 , (290) 

a 3 (s, Sl ) = <%(s,8i)+a$ (s,si), (291) 

[3-3s-3s 2 +4si(sx + l)] 2 



c#(s, Sl ) = 



2(2s + 3) 2 
(s - 1)[3 + 3s - 3s 2 + 4si(si + l)] 2 



of (s,si) 



2s(2s + l) 2 



[s(s + 2)-4s 1 (s 1 + l)] 
2(s + l)(s + 3)(2s + l) 2 (2s + 3) 2 
„ Q3!o(g)+4si(si + l)a^ 1 (s) 
(2s + 5)(3s + l) 



(292) 



a* (s) = (s+l)(s + 3)(51 + 114s 



a 4 (s, s x ) 



06 ( s > s 



Oe (s,Sl) 



+209s 2 + 302s 3 + 164s 4 + 24s 5 ), (293) 



afi(s) = 129 + 318s + 395s 2 



+442s 3 + 300s 4 + 72s 5 , 
(4f - 3) 
8s ' 

3[a 5 ,o(s) + 4si(si + 1)] 



as(s ' Sl) ~" 4s(2s + l)(2s + 3) ' 

O5.0 («) = 3 + 3s - 5s 2 - 4s 3 , 
06(s,s x ) = a^(s,si) +a%(s,si), 

[3-3s-3s 2 +4si(s! + l)] 2 



(294) 

(295) 

(296) 
(297) 



8(2s- l)(2s + 3) 2 
(s - 1)[3 + 3s - 3s 2 + 4si(si + l)] 2 



8s(2s- 3)(2s+ l) 2 
1 



(298) 



4(2s - 3) (2s - l)(2s + l) 2 (2s + 3) 2 

En = 0<nW[4ai(«l + l)]" 

(2s + 5)(4s + l)(4s + 3) ' 
Oeo( s ) = s(-216 + 837s + 5052s 2 + 3521s 3 



T2414s 4 - 21876s 5 - 7464s 6 



+6720s 7 + 5376s 8 + 1024s 9 ) , (299) 
a^s) = 3(81 + 102s + 10s 2 

+ 196s 3 + 296s 4 + 80s 5 ), (300) 



(301) 



a 6,2( s ) = 189 + 258s -100s 2 + 152s 3 
+640s 4 + 256s 5 , 
_ [s(s + 2)-4s 1 (s 1 + l)] 
n ' 11 4(2s+l) 3 (2s + 3) 3 (2s + 5) 

x U Sl (sx + 1)(-1 + 38s + 60s 2 + 24s 3 ) 

+ (s + 1)(3 + 67s + 94s 2 + 44s 3 + 8s 4 ) 

(302) 

By expanding the solutions in the crystal representa- 
tion to second order in the Jj , we have explicitly checked 
these expressions for s% = 1/2, s = 1, and for s± = 1, 
s = 1,2. We note that for si = 1/2, only a\ and 04 are 
non- vanishing. 

Two-level thermodynamic coefficients ai, s 

Here we present the expression for oi )S to 0(Jj/J) 2 
appearing in Eqs. and (j4TTl -(l4~4l in the text. We 

have 



Oi 



J2b n ( s , Sl )fi 2 \e, ( f>) + o(j j /j) 3 



(303) 
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where the /„(#, 0) are given by Eqs. H216(l - I|221|) . ancl 
(2s - l) 2 



h(s,si) 

M s ; s i) 

h(s, si) 



2s 

(2* -I) 2 

(2s — l) 2 
2s 



6(s-l), 
a SjSl 6(s - 1) 



(304) 
(305) 



6 4 (s,si) 

&5(s,Sl) 



< Sl 6( S -l) + 6^( SlSl ), (306) 

s(27 + 44s + 27s 2 + 6s 3 ) 
~ 6(s + l) 2 (s + 3) 2 (2s + 3) 2 (2s + 5) 

x([4si(si + l)-s(s + 2)] 

x[4s 1 (s 1 + l)-(s + l)(s + 3)]), (307) 
(2s -1) 



8s 
(2s - 1 
IT" 



-e(s-i), 

■a StS1 Q(s - 1), 



(308) 
(309) 



be(s, Sl ) = i^-ila^e^-^ + fe^s,^), (310) 



x[s(s-2)-4 Sl (si + l)] 

s6 64 (s)[(s(s + 2)-4s 1 (s 1 + l)] 
72(2s + l)(2s + 3) 2 (2s + 5)(4s + 3) 2 
x[(s + l)(s + 3)-4s 1 (s 1 + l)], (311) 
6 M ( S ) = 369 + 1011s + 1420s 2 + 1076s 3 



+416s 4 + 64s 5 , 



where 8(x) = 1 for x > 0, <d(x) = for x < is the 
Heaviside step function. 



[1] 
[2] 
[3] 
[4] 

[5] 



Electronic address: efremov@theory.phy.tu-dresden.de 
Electronic address: klemm@phys.ksu.edu 
R. Sessoli, D. Gatteschi, A. Caneschi, and M. A. Novak, 
Nature (London) 365, 141 (1993). 

J. R. Friedman, M. P. Sarachik, J. Tejada, and R. Ziolo, 
Phys. Rev. Lett. 76, 3830 (1996). 

M. N. Leuenberger and D. Loss, Nature (London) 410, 
789 (2001). 

W. Wernsdorfer, T. Ohm, C. Sangregorio, R. Sessoli, D. 
Mailly, and C. Paulsen, Phys. Rev. Lett. 82, 3903 (1999). 
W. Wernsdorfer and R. Sessoli, Science 284, 133 (1999). 



[6 
[T 

[8] 

[9] 

[10 

[11 

[12 
[13 

[14 

[is; 

[16 

[17; 

[18 



[19 
[20 



(312) [21 



[22 

[23 

[24 

[25 
[26 

[27 

[28 

[29 



D. Zipse, J. M. North, N. S. Dalai, S. Hill, and R. S. 
Edwards, Phys. Rev. B 68, 184408 (2003). 
M. Affronte, A. Cornia, A. Lascialfari, F. Borsa, D. Gat- 
teschi, J. Hinderer, M. Horvatic, A. G. M. Jansen, and 
M.-H. Julien, Phys. Rev. Lett. 88, 167201 (2002). 
O. Waldmann, J. Schiilein, R. Koch, P. Miiller, I. Bernt, 
R. W. Saalfrank, H. P. Andres, H. U. Gudel, and P. Al- 
lenspach, Inorg. Chem. 38, 5879 (1999). 
O. Waldmann, R. Koch, S. Schromm, J. Schiilein, P. 
Miiller, I. Bernt, R. W. Saalfrank, F. Hempel, and E. 
Balthes, Inorg. Chem. 40, 2986 (2001). 
H. Nakano and S. Miyashita, J. Phys. Chem. Solids 63, 
1519 (2002). 

D. A. Tennant, S. E. Nagler, A. W. Garrett, T. Barnes, 

and C. C. Torardi, Phys. Rev. Lett. 78, 4998 (1997). 

H. U. Giidel, Neutron News 7, 24 (1996) 

A. W. Garrett, S. E. Nagler, D. A. Tennant, B. C. Sales, 

and T. Barnes, Phys. Rev. Lett. 79, 745 (1997). 

D. V. Efremov and R. A. Klemm, Phys. Rev. B 66, 

174427 (2002). 

F. Le Gall, F. Fabrizi de Biani, A. Caneschi, P. Cinelli, A. 
Cornia, A. C. Fabretti, and D. Gatteschi, Inorg. Chim. 
Acta 262, 123 (1997). 

Y. Shapira, M. T. Liu, S. Foner, R. J. Howard, and W. 

H. Armstrong, Phys. Rev. B 63, 094422 (2001). 

Y. Shapira, M. T. Liu, S. Foner, C. E. Dube, and P. J. 

Bonitatebus, Jr., Phys. Rev. B 59, 1046 (1999). 

K. L. Taft, C. D. Delfs, G. C. Papaefthymiou, S. Foner, 

D. Gatteschi, and S. J. Lippard, J. Am. Chem. Soc. 116, 

823 (1994). 

J. D. Walker and R. Poli, Inorg. Chem. 29, 756 (1990). 
J. A. Bertrand, J. L. Breece, and P. G. Eller, Inorg. 
Chem. 13, 125 (1974). 

R. Tiron, W. Wernsdorfer, D. Foguet-Albiol, N. Aliaga- 
Alcalde, and G. Christou, Phys. Rev. Lett. 91, 227203 
(2003). 

J. M. North, N. S. Dalai, D. Foguet-Albiol, A. Vinslava, 
and G. Christou, Phys. Rev. B 69, 174419 (2004). 
O. Waldmann, J. Hassmann, P. Miiller, D. Volkmer, U. S. 
Schubert, and J.-M. Lehn, Phys. Rev. B 58, 3277 (1998). 
J. J. Borras-Almenar, J. M. Clemente-Juan, E. coronado, 
and B. S. Tsukerblat, Inorg. Che m. 38, 6081 (1999). 
D. V. Efremov and R. A. Klemm, jcond-ma t /0409168 
R. M. White, Quantum Theory of Magnetism (McGraw- 
Hill, New York, 1970), Chapters 1, 5, 8. 
K. Park, M. A. Novotny, N. S. Dalai, S. Hill, and P. A. 
Rikvold, Phys. Rev. B 66, 144409 (2002). 
R. A. Klemm and J. R. Clem, Phys. Rev. B 21, 1868 
(1980); R. A. Klemm, SIAM J. Appl. Math. 55, 986 
(1995). 

H. Goldstein, Classical Mechanics, (Addison- Wesley, 
Reading, MA, 1965), p. 109. 



